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INTRODUCTION 

The  work  summarized  in  Ref.  1 has  shown  that  a functional  poly- 
nomial input-output  model  has  considerable  practical  promise  as  a 
unifying  concept  for  the  interpretation  of  ship  added  resistance  and 
possibly  related  problems.  In  particular,  in  this  and  previous  work  it 
has  been  demonstrated  that;  a)  both  a linear  and  a quadratic  frequency 
response  function  can  be  derived  from  experiments  in  irregular  waves,  as 
well  as  from  experiments  in  regular  waves;  b)  synthesis  of  the  mean  added 
resistance  and  of  resistance  spectra  can  be  carried  out  by  use  of  these 
frequency  response  functions,  and;  c)  time  histories  of  added  resistance 
can  be  synthesized  for  an  irregular  wave  input  through  use  of  the  time 
domain  representation  of  the  linear  and  quadratic  level  non-linear  fre- 
quency response  functions. 

All  of  the  just  cited  work  was  empirical.  The  particular  input- 
output  model  was  hypothesized  at  the  outset  and  the  implications  of  the 
mode)  were  checked  as  far  as  possible  with  experiment.  The  idea  that 
there  exist  linear  and  quadratic  frequency  response  functions  followed 
from  the  (non-physical)  model  assumed.  The  input-output  model  provides 
no  detailed  physical  basis.  In  the  context  of  added  resistance,  a linear 
frequency  response  function  may  be  thought  of  as  the  relation  between 
oscillatory  surge  exciting  force  and  wave  elevation.  There  are  a number 
of  hydromechanical  approaches  to  the  estimation  of  this  function  so  that 
the  concept  could  be  accepted  as  having  an  identifiable  physical  base. 

On  the  other  hand,  methods  of  hydromechanical  computation  for  the  entire 
quadratic  response  function  were  lacking.  Only  a very  special  portion 
of  the  quadratic  frequency  response  function  could  be  estimated  from 
hydromechanical  cons iderat ions . This  problem  was  addressed  in  the  work 


1.  Dalzell,  J.F.,  "Application  of  the  Functional  Polynomial  Model  to 
the  Ship  Added  Resistance  Problem,"  Eleventh  Symposium  on  Naval 
Hydrodynamics,  London,  March  1976. 
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of  Ref.  2 , and  in  this  latter  work  the  relationship  between  the  hydro- 
mechanics of  added  resistance  and  the  "quadratic  frequency  response 
function"  of  the  input-output  model  was  clarified.  Reasonably  good 
qualitative,  and  fair  quantitative  agreement  was  achieved  between  results 
of  analysis  and  experiment.  Thus  the  work  of  Ref.  2,  by  providing  a 
physical  model,  lends  increased  confidence  to  the  use  of  the  two-term 
functional  polynomial  model  in  analysis  and  interpretat ion  of  experiment. 

Insofar  as  ship  model  dynamics  in  waves  is  concerned,  the  general 
experimental  problem  is  to  identify  the  pertinent  parameters  or  functions 
given  a wave  input  and  the  observed  response.  An  input-output  model  of 
some  sort  is  ordinarily  required  if  an  endless  series  of  experiments  is 
to  be  avoided. 

Historically,  the  advent  of  the  linear  input-output  model  for  some 
ship  dynamics  problems  had  the  effect  of  expanding  the  possibilities  in 
experimental  ship  dynamics  work.  Once  it  was  shown  that  the  linear  model 
was  a reasonable  engineering  approach  for  the  prediction  of  ship  response 
under  realistic  (random)  conditions,  three  experimental  techniques  were 
admissible.  Regular  wave  experiments  took  on  a slightly  different  mean- 
ing, the  option  became  available  for  the  interpretation  of  experiments 

in  irregular  waves  by  spectral  analysis,  and  the  development  of  transient 
3* 

test  techniques  followed  shortly  thereafter.  Of  the  three  techniques 

th*  first  two  are  the  most  widely  used,  although  the  transient  technique 

4 5* 

has  undergone  additional  development  in  recent  years,  Takezawa,  et  al  ’ 


■sV 

2.  Dalzell,  J.F.  and  Kim,  C.H.,  "Analytical  Investigation  of  the  Quadratic 
Frequency  Response  for  Added  Resistance,"  S IT-DL-76- 1 878,  Davidson 
Laboratory,  Stevens  Institute  of  Technology,  August  1976. 

3.  Davis,  M.C.  and  Zarnick,  E.E.,  "Testing  Ship  Models  in  Transient 
Waves,"  Fifth  Symposium  on  Naval  Hydrodynamics,  Bergen,  1964. 

4.  Takezawa,  S.  and  Takekawa,  M.,  "Advanced  Experimental  Techniques  for 

Testing  Ship  Models  in  Transient  Water  Waves:  Part  I,  The  Transient 

Test  Technique  on  Ship  Motions  in  Waves,"  Eleventh  Symposium  on 
Naval  Hydrodynamics,  London,  1976. 

5.  Takezawa,  S.  and  Hirayama,  T.,  "Advanced  Experimental  Techniques  for 

Testing  Ship  Models  in  Transient  Water  Waves:  Part  II,  The  Controlled 

Transient  Water  Waves  for  Using  in  Ship  Motion  Tests,"  Eleventh 
Symposium  on  Naval  Hydrodynamics,  London,  1976. 
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Given  this  history  and  the  evidence  thus  far  advanced  that  the 
two  term  functional  polynomial  is  a realistic  input-output  model  for 
added  resistance,  it  is  natural  to  inquire  into  the  extent  to  which  tech- 
niques paralleling  the  three  available  for  the  linear  case  are  available 
for  routine  experimental  determination  of  the  quadratic  frequency  response 
funct ion. 

A regular  wave  technique  has  been  used  for  some  time  in  the  determ- 
ination of  mean  added  resistance  operators.  In  this  technique  the  total 
mean  resistance  of  a model  in  regular  waves  is  measured  and  the  added 
resistance  operator  is  derived  from  the  difference  between  this  observation 
and  the  model  resistance  in  calm  water.  Multiple  runs  are  required  to 
cover  a range  of  wave  frequency  and  it  has  been  found  by  some  that  frequent 
re-running  of  calm  water  tests  is  advisable.  Many  of  the  basic  problems 
with  the  technique  appear  to  be  related  to  experimental  accuracy.  The  mean 
added  resistance  is  usually  small  in  comparison  with  calm  water  resistance, 
and  both  are  often  very  small  in  comparison  with  the  time  dependent  com- 
ponent . 

Because  the  mean  added  resistance  operator  corresponds  only  to  a 

special  portion  of  the  quadratic  frequency  response  function,  the  existing 

regular  wave  technique  is  not  sufficient  if  experimental  estimates  of  the 

entire  function  are  required.  The  principle  of  the  required  regular  wave 

technique  has  been  indicated  in  Reference  1.  However,  this  technique  has 

not  so  far  been  used  in  practice.  In  brief,  to  produce  results  defining 

the  entire  quadratic  response  function  for  added  res i s tance, an  experimental 

technique  involving  dual  harmonic  excitation  is  required;  that  is,  two 

superimposed  regular  waves.  This  technique  corresponds  to  the  basic 

interpretation  of  the  quadratic  frequency  response  function  given  in  Ref.  1. 

It  would  involve  selection  of  pairs  of  regular  wave  frequencies  such  that 

their  superposition  in  the  encounter  domain  is  periodic,  as  well  as  the 

subsequent  harmonic  analysis  of  the  added  resistance  response.  The  general 

nature  of  the  quadratic  frequency  response  function  and  the  principle  of 

the  regular  wave  technique  imply  that  if  N experimental  runs  in  regular 

waves  are  adequate  to  describe  the  important  part  of  the  added  resistance 

2 

operator,  then  about  N / 2 runs  might  be  required  for  a reasonable  definition 
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of  the  entire  function.  Thus  a 10  run  per  model  speed  mean  added  resist- 
ance test  program  might  turn  into  a 50  run  per  speed  program  for  the 
general  quadratic  response. 

With  respect  to  the  possibilities  of  identifying  the  linear  and 
the  quadratic  frequency  response  functions  from  data  obtained  in  irreg- 
ular waves,  the  techniques  are  more  or  less  available.  Ordinary  spectral 
and  cross-spectral  analyses  suffice  for  the  linear  function,  and  a tech- 
nique called  cross-bi-spectral  analysis  has  been  developed  for  the 
quadratic  function  ( Ref . 1 ) . In  principle  the  cross-bi-spectral  analysis 
technique  is  quite  attractive  in  the  sense  of  increasing  efficiency  over 
that  for  regular  wave  experiments.  The  technique  is  based  on  a correla- 
tion with  the  wave  fluctuations  of  the  fluctuations  in  resistance.  Thus 
this  technique  is  not  sensitive  to  the  value  of  calm  water  resistance,  and 
in  fact  calm  water  resistance  need  not  be  known  at  all.  Again  in  principle, 
either  the  mean  added  resistance  or  the  entire  quadratic  frequency  response 
may  be  derived  from  the  same  sample  of  irregular  wave  input  and  response, 
the  difference  being  the  amount  of  computer  time  brought  to  bear. 

In  practice  however,  the  technique  has  a quite  serious  deficiency.  It  is 
that  for  reasonable  precision  of  estimates  an  order  of  magnitude  larger 
sample  is  presently  required  than  is  customary  in  the  spectral  identifica- 
tion of  the  linear  response.  In  many  cases  it  is  not  possible  to  accumulate 
sufficient  sample  for  the  linear  identification  process  in  one  pass  up  a 
towing  tank.  Thus  it  appears  that  the  cross-bi-spectral  analysis  technique 
will  always  require  as  many  or  more  test  runs  as  the  existing  regular  wave 
technique  if  only  the  mean  added  resistance  operator  is  required  — although 
fewer  for  a determination  of  the  entire  function.  The  mean  added  resistance 
operator  is  presently  the  only  part  of  the  quadratic  frequency  response 
function  which  finds  routine  use,  and  it  is  to  be  expected  that  the  situa- 
tion will  not  soon  materially  change.  Accordingly,  it  does  not  now  appear 
that  the  technique  involving  cross-bi-spectral  analysis  of  irregular  wave 
test  data  will  offer  savings  in  facility  time  over  that  now  required  for 
routine  added  resistance  experiments. 

In  parallel  with  the  possibilities  for  the  linear  case,  the  third 
experimental  technique  involves  the  analysis  of  transient  wave  pulses. 


I 

I 
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Conceptually,  the  transient  wave  pulse  would  involve  all  wave  frequencies 
of  interest,  and  the  wave  elevation  would  be  zero  or  nearly  so  at  begin- 
ning and  end  of  experiment.  Thus  the  observed  resistance  transient, 
while  containing  response  at  all  frequencies,  should  equal  calm  water 
resistance  at  beginning  and  end  of  the  experiment.  The  ideal  advantages 
are  thus  similar  to  those  for  the  cross-bi-spectral  analysis  technique  in 
that  all  the  required  information  is  embedded  somehow  in  a single  observa- 
tion. The  difficulty  with  the  idea  is  essentially  that  the  existing 
theory  for  the  quadratic  input-output  model  does  not  contain  a clear 
indication  of  how  to  proceed  with  the  analysis. 

It  was  accordingly  the  objective  of  the  present  work  to  investi- 
gate the  feasibility  of  a wave  pulse  technique  for  added  resistance.  For 
the  practical  reasons  cited  previously  the  emphasis  was  to  be  upon  ident- 
ification of  the  mean  added  resistance  operator  rather  than  the  entire 
quadratic  frequency  response  function,  and  the  methods  to  be  employed  were 
to  involve  digital  computer  rather  than  physical  experiments. 
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THE  TRANSIENT  EXPERIMENT 

In  order  to  fix  some  terminology,  various  stages  of  the  hypotheti- 
cal transient  experiment  for  head  seas  are  outlined  in  Figure  1.  In  the 

figure  the  condition  of  tank  and  model  are  indicated  for  five  significant 

stages  of  one  experimental  run.  Time  (t)  will  be  considered  to  be  posi- 
tive after  the  wave  generator  is  started.  Thus  there  is  shown  at  the  top 
of  the  figure  for  t < 0 a hypothetical  towing  tank  which  has  a wave  gen- 
erator at  one  end,  an  absorber  at  the  opposite  end,  and  a model  sitting 
still  in  preparation  for  the  run.  At  this  stage  the  water  is  assumed  to 

be  quiet  everywhere.  Position  in  the  tank  is  denoted  by  X,  which  is  zero 

at  the  wavemaker  and  positive  in  the  direction  of  wave  propagation. 

At  time  equal  to  zero  the  wave  generator  begins  to  generate  a long- 

crested  transient  wave,  and  typically  at  some  subsequent  time  the  model 

acceleration  begins  so  that  at  the  second  stage  (t=ts)  the  model  is 

proceeding  with  constant  velocity,  U,  at  tank  position  X=Xs.  The  model 

is  assumed  to  be  in  still  water  at  this  stage  while  the  wave  transient 

at  the  other  end  of  the  tank  is  propagating  toward  it.  In  the  third 

stage  (t=tst)  the  model  enters  the  wave  pulse  at  tank  position  X=X5t* 

At  the  fourth  stage  (t=t  ) the  model  leaves  the  pulse  at  position  X , 

et  et 

and  continues  at  constant  speed  in  calm  water  until  time  t£  (position  Xg) 

when  deceleration  begins.  The  useful  part  of  the  experiment  is  between 

times  t$  and  tg,  during  which  an  observer  on  the  model  sees  calm,  then 

rough,  and  then  finally  calm  water  at  the  end.  Because  of  the  assumption 

of  constant  speed,  the  position  of  the  model  in  the  tank  is  known  during 

this  time  from  a knowledge  of  X , t and  U. 

a s'  s 

Several  of  the  tacit  assumptions  in  the  above  deserve  cormient. 

One  is  that  of  a long-crested  wave  pulse  or  transient.  The  input-output 
model  has  not  been  extended  to  the  short-crested  case.  For  the  purposes 
of  the  present  study  this  is  not  a serious  restriction  since  none  of  the 
other  techniques  described  in  the  introduction  have  been  adapted  for  this 
case,  and  in  fact  there  have  as  of  the  present  been  very  few  experiments 
of  any  type  in  short-crested  seas. 
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I 

Another  comment  which  can  be  made  is  about  the  subject  of  heading. 

Clearly,  if  there  is  some  component  of  wave  group  velocity  in  the  direc- 
tion of  model  travel  there  may  be  an  ambiguity  which  would  require  very 
special  handling,  or  the  model  may  never  pass  out  of  the  transient  at 
all  — at  least  not  within  the  confines  of  a towing  tank.  For  example, 
in  following  seas  the  model  might  be  started  near  the  wave  maker  before 
the  pulse  is  generated.  If  all  wave  components  travel  at  group  velocity 
greater  than  model  velocity,  the  wave  pulse  will  overtake  and  pass  the 
model.  Alternately  if  the  model  speed  is  high  enough  it  might  be  possible 
for  the  model  to  overtake  and  pass  the  pulse.  In  between  these  extremes 
it  is  extremely  difficult  to  envision  a reasonable  following  seas 
transient  experiment  when  the  model  speed  and  the  mean  group  velocity  of 
the  wave  components  are  about  the  same.  For  purposes  of  the  present  study 
the  problems  about  all  model  headings  excepting  head  seas  were  ignored  on 
the  grounds  that  the  feasibility  of  even  the  head  sea  test  was  not  clear, 
and  that  the  head  sea  case  is  usually  considered  the  most  important  for 
addt  r istance. 

The  assumption  of  constant  model  velocity  implies  that  the  model 
is  restrained  in  surge  and  that  resistance  is  measured  through  some  sort 

of  force  balance.  This  is  a system  not  always  used  in  practice.  It  was 

6* 

retained  in  the  present  case  because;  a)  it  is  feasible  to  do  b)  it  has 
7 8* 

been  found'*  that  the  mean  added  resistance  for  constant  speed  is  not 
significantly  influenced  by  surge  restraint;  and  c)  it  seemed  probable 
that  a 'technique  developed  for  this  case  might  be  adaptable  to  other  con- 
ditions of  surge  restraint. 

V? 

6.  Dalzell,  J.F.,  "Application  of  Cross-B i -Spectra  1 Analysis  to  Ship 
Resistance  in  Waves,"  SIT-DL-72-l606,  AD  7^+9 102,  Davidson  Laboratory, 

Stevens  Institute  of  Technology,  May  1972. 

7.  Sibul,  O.J.,  "Constant  Thrust  vs.  Constant  Velocity  Method  for 
Resistance  Measurement  in  Waves,"  University  of  California,  Berkeley, 

Report  NA-71-1,  June  1971. 

8.  Journee,  J.M.J.,  "Motions,  Resistance  and  Propulsion  of  a Ship  in 
Longitudinal  Regular  Waves,"  Report  428,  Laboratorium  voor 
Sheepshydromechanica,  Technische  Hogeschool  Delft,  May  1976. 
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WAVE  PULSES 


Theory 

It  was  the  basic  objective  of  the  work  to  examine  feasibility  of 
techniques  which  might  be  employed  in  towing  tanks,  and  to  check  any 
promising  scheme  with  a computer  "experiment."  Accordingly  it  was  neces- 
sary to  begin  with  a consideration  of  the  type  of  wave  transient  which  is 
physically  possible,  and  to  generate  a few  examples  for  later  use. 
References  3 and  5 indicate  clearly  that  because  of  the  dispersive  proper- 
ties of  waves,  it  is  not  possible  to  produce  anything  closely  resembling 
the  ideal  impulse,  or  the  isolated  step  function,  etc.,  which  are  basic 
to  much  of  the  theory  of  transient  response  as  shown  in  the  textbooks  of 
control  and  electronics. 

The  linear  theory  of  the  type  of  wave  pulse  which  is  possible  is 
well  developed,  is  treated  in  References  3 and  5>  and  will  be  used  herein. 
It  is  apparent  from  these  references  that  the  linear  theory  and  experi- 
ment correlate  well  so  long  as  average  local  wave  steepnesses  are  not 
such  as  to  produce  extensive  breaking. 

Following  the  developments  in  References  3 and  5,  the  transient 
wave  elevation  at  position  X in  the  tank  will  be  denoted  h^(t) . It  is 
assumed  that  h^(t)  is  absolutely  integrable,  that  is: 

J I hx(t)  | dt  < co 

where  the  convention  is  followed  that  omission  of  limits  signifies  limits 
of  plus  and  minus  infinity.  (This  convention  will  be  followed  throughout 
the  present  report.)  On  the  basis  of  this  assumption,  the  transient  has 
a Fourier  transform  or  complex  spectrum,  (uj)  , and  the  transform  pair 
relating  the  transient  and  its  spectrum  may  be  defined  as: 
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hx(t)  = -57  JHx(uj)  Exp[  iujt]duj  (1) 

= j\(t)Exp[-iaJt]dt  (2) 

where  uj  denotes  circular  wave  frequency. 

In  the  practical  production  of  a wave  transient,  the  wavemaker  is 
given  a specified  sequence  of  control  signals  so  as  to  produce  a desired 
transient  wave  elevation.  For  present  purposes  the  dynamics  of  the 
wavemaker  may  be  neglected  and  it  will  be  assumed  that  the  transient  wave 
at  the  wavemaker  is  specified.  Thus  the  wave  at  X=0  (Figure  1)  will  be 
taken  to  be: 

hoW  = 57  iH0^— ^ Expt  iujfJduj  (3) 

where 

Hq(u))  is  the  complex  spectrum  of  the  wave  pulse. 

There  are  two  main  constraints  on  hQ(t) . The  first  is  that  it  be  zero 
for  t < 0;  that  is,  hQ(t)  must  be  physically  realizable  according  to 
the  time  conventions  shown  in  Figure  1.  The  second  is  that  the  local 
wave  steepnesses  be  "small." 

Following  the  linear  theory  of  wave  propagation^'"’,  the  frequency 
response  function  connecting  the  wave  elevation  at  X=0  (Eq.3)  and  wave 
elevation  elsewhere  in  the  tank  is  given  by: 

Exp[-i  | u)  | u)X/g] 

which  simply  means  that  each  wave  component  is  assumed  to  propagate  with- 
out change  in  amplitude.  From  this,  the  spectrum  of  the  wave  transient 
at  position  X becomes: 

Hx(uj)  = HQ((u)Exp[-i  | u>  | #/g]  (4) 

and  the  transient  wave  itself  at  position  X is: 

hx(t)  JHo(uj)  Exp[iujt-i  | a)  | UJX/g]du|  (5) 
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Also  following  from  the  linear  assumption  is  that  the  integral  square  of 
the  transient  is  invariant  with  tank  position: 

Ex  = Jhx^  dt  " ^ J I f d«i  * constant  (6) 

in  which  the  last  relation  is  obtained  by  application  of  the  Parseval 
Theorem. 

The  invariance  of  the  integral  square  suggests  a useful  change  in 
notation  for  the  complex  spectrum  of  the  transient  at  the  wavemaker: 

Ho(m)  = A(uj)Exp[icp(o))]  (7) 

in  which  A(u>)  and  cp(u>)  are  real 
and : 

A(-uj)  = A (iu) 
cp(-co)  = -cp(a;) 

Then  the  wave  transient  at  position  X becomes: 

h (t)  = JA(uj)  Exp[  iu)t+i9(u>)  - i | u)  | ujX/g]doi  (8) 

In  the  transient  techniques  described  in  References  3,4,5  the  amplitude 
part  of  the  spectrum  of  the  wave  at  the  wavemaker,  A(uj),  is  specified 
so  as  to  produce  wave  excitation  in  all  the  frequencies  of  interest. 

The  phase  part,  T>(<u) , is  used  to  control  the  wave  pulse  in  the  run  area. 

It  is  shown  in  Reference  5 that  the  wave  pulse  of  shortest  dura- 
tion has  a constant  phase  spectrum.  Accordingly,  for  the  shortest 
duration  of  pulse,  it  is  necessary  to  specify  a part  of  <P(uj)  so  that 
the  | u)  | U)X/g  term  in  Eq.  (8)  is  cancelled  out  at  the  position  of  interest. 
Thus  the  phase  function  for  the  transient  at  the  wavemaker  may  be  speci- 
fied as  follows: 

cp(uj)  = Cn  sgn(u))+  | uj  | jf^/g  - u)tm  (9) 

In  Eq.  (9)  C is  an  arbitrary  constant.  The  second  term  is  a phase  lead 
which  will  cancel  the  similar  phase  lag  in  Eq.  (8)  at  position  X^.  The 
last  term  is  a phase  lag  which  corresponds  to  a constant  time  delay  in 
the  time  domain.  It  must  be  included  in  the  phase  prescription  as  a way 
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to  insure  that  the  transient  at  the  wavemaker  will  be  physically  realiz- 
able. What  is  accomplished  by  the  above  assumption  is  to  organize  the 
transient  wave  at  the  wavemaker  in  such  a way  that  about  an  equal  time  is 

necessary  for  all  component  waves  to  reach  position  X in  the  tank. 

m 

At  positions  closer  to  the  wavemaker  the  transient  pulse  is  converging, 
and  at  positions  further  away  than  Xm  the  transient  is  diverging. 

In  practical  testing  in  the  usual  towing  tank  a desirable  mid- 
point of  the  data  taking  portion  of  the  run  is  usually  known  from  the 
acceleration  and  deceleration  characteristics  of  the  carriage  drive,  and 
from  the  positions  of  the  inevitable  obstructions  at  either  end  of  che 
tank.  Thus  the  tank  position,  X^,  in  Eq.  (9)  might  reasonably  be  speci- 
fied as  the  mid-point  of  the  data  run,  and  would  ordinarily  be  determined 
by  other  than  wave  making  considerations.  In  the  terminology  of  Figure  1: 

Xm  ’ (Xs  * Xe)/2  <l0> 

Thus  substituting  Eq.  (9)  into  Eq.  (8)  there  results  a further 
expression  for  the  wave  pulse  at  any  tank  position  X: 

hx(t)  = jA(ui)  Exp[  iu|  (t-tm)  + iCrrsgn(uj)-i  ] u;  | uj(X -Xm) /g]dw 

= - J A(u))cos[uj(t-tm)+CTT  - U)1 2  (X-Xm)/g]du) 
o 

oo 

= - f A(u))cos[CTT-u)t  - a)2  (X-X  )/glcosu)t  duj 
tt  « - — m — ' m 3 — — 

o 

1 00 

- — J A(uj)sin[Cn-u)tm  - u;2  (X-X  )/g]sinu;t  duj  (II) 

o 

In  the  present  case  the  encountered  wave  transient,  h (t),  is 
desired.  It  may  be  assumed  that  timing  of  model  acceleration  relative  to 
wavemaker  start  can  be  made  such  that  the  model  will  arrive  at  position  X 

m 

at  time  (tm-6) . For  constant  model  velocity,  U,  the  position  of  the  model 
in  the  tank  may  then  be  written  as: 

X ■ Xm  - U'*-Vf)  02> 

In  the  practical  case  Eq.  (12)  is  valid  only  within  the  range  of  time 
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(ts  ^ t ^ t ).  Now  substituting  Eq.  (12)  into  Eq.  (11)  there  results  an 
expression  for  the  encountered  transient: 

hg(T+tm)  = -^  ^(u))  Exp[  ix(oj  + | a;  | ujU/g)  ]duj  (13) 

where:  T=t-t 

m 

B(uj)  = Br(uj)  + i B j ( uj) 

with: 

Br(uj)  = A(uj)  cos[Crrsgn(uj)  + | uj  | ujU6/g] 

B | (u>)  = A(uj)sin[Cnsgn(uj)  + | u)  | ujU6/g]  (14) 

With  respect  to  the  practical  experiment,  Eq . (13)  is  valid  only  in  the 
range  (ts“tfn)  ^ t s (t  -tm)  . The  time  variable  change  corresponds  to  the 
normal  practice  in  experiments  of  counting  time  with  reference  to  some 
arbitrary  point  in  the  constant  speed  portion  of  the  run.  The  function 
B(uj)  contains  only  parameters  which  are  prescribed  by  programming  of  the 
wave  machine,  and  the  timing  of  model  start.  The  function  of  uj  in  the 
exponential  of  Eq.  (13)  may  be  recognized  as  the  encounter  frequency  for 
head  seas,  uj.  The  expression,  Eq . (13),  may  thus  be  transformed  into  an 
integration  over  U).  As  such  it  is  in  the  form  of  a Fourier  Transform. 

So  long  as  matters  are  arranged  so  that  h (t)  is  zero  outside  the  range  of 
practical  validity  of  Eq.  (13),  the  inverse  transform  of  the  encountered 
transient  will  be  valid  in  the  sense  that  there  will  be  no  difference 
between  the  practical  experiment  and  the  hypothetical  situation  in  which 
the  model  proceeds  at  constant  velocity  for  an  infinite  range  of  x. 

Applying  the  head  sea  transformation: 

uj  = uj  + ! uj  | u)U/g 

and  defining: 

K(oj)  = fjj  (-1  + /1+Wg) 

J(uj)  = l//l+4ujU/g  (15) 
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00 

h (T+f  ) = £ J*  B^(K(oi) ) J (to)  cos  (ujt)  dm 
o 

1 00 

"■  - — J*  B|  (K(uj))j(u))sin(a)T)du)  (16) 

o 

The  effective  Fourier  spectrum  of  the  encountered  transient  is  B (K(m)  ) J (u))  . 

From  Eq.  (14)  : 

| B(K(iu)  ) J(uj)  | = A(K(m))J(u>)  (17) 

so  that,  as  would  be  expected,  the  amplitude  part  of  the  Fourier  spectrum 
is  distorted  by  the  transformation. 

Example  Evaluations 

As  noted,  some  examples  of  physically  possible  wave  pulses  were 
desired.  The  easiest  way  to  insure  realism  was  to  evaluate  Eq.  (11)  with 
parameters  appropriate  to  a given  tank  and  model  size.  The  Davidson 
Laboratory  Tank  No.  3 (300lxl2'x6l)  was  selected  as  the  hypothetical  tank, 
and  the  hypothetical  experiment  was  assumed  to  involve  a five  foot  ship 
model  at  a speed  corresponding  to  Froude  Number  0.15.  ( 

The  first  function  to  specify  is  A(uj)  . It  was  assumed  on  the  basis 
of  prior  experimental  results^  that  wave  lengths  between  2.5  and  15  feet  ' 

(1/2  to  3 times  model  length)  would  include  the  important  parts  of  the 
resistance  response.  It  was  further  assumed  that  for  the  identification 
purposes  envisioned  it  would  be  of  advantage  to  make  the  encountered  spec- 
trum have  a roughly  constant  modulus  ( Eq .17)  in  the  encounter  frequency 
range  which  corresponds  to  the  wave  length  range  just  cited.  The  first 
step  in  constructing  A (u»)  was  thus  to  make  A (uj)  = A ( K ( U)) ) = P/J(tu)  where  P 
is  a constant  which  controls  the  magnitude  of  the  generated  transient. 

The  model  speed  assumed  in  computing  encounter  frequencies  corresponded  to 
a Froude  Number  of  0.15  for  a five  foot  model.  The  resulting  form  for  A(uj) 
is  shown  in  Figure  2.  The  wave  frequency  range  of  interest  was  roughly 
3.5  to  9 radians/second.  In  this  range  the  function  is  as  noted  above. 

At  each  end  of  this  range  a smooth  transition  is  made  to  zero  over  a 
one  radian/second  frequency  range. 
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FIGURE  2 FORM  OF  SPECIFIED  MODULUS 
OF  SPECTRUM  OF  WAVEMAKER 
TRANSIENT  FOR  FROUDE  NUMBER  0.15 
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The  next  parameters  required  by  Eq.  (ll)  are  Xm  and  t . In  D.L. 

Tank  3 the  mid-point  of  the  data  taking  area  is  usually  150  feet  from  the 
wavemaker,  and  this  distance  was  chosen  for  X^,  the  designed  concentra- 
tion point.  The  time  delay  was  chosen  to  be  a bit  more  than  the  time 
required  for  the  wave  component  of  highest  frequency  to  travel  150  feet 
at  its  group  velocity.  The  number  used  was  t = 105  seconds. 

The  constant,  C,  in  Eq.  (11)  controls  the  symmetry  of  the  wave 
transient  at  the  point  of  maximum  concentration.  If  the  constant  is  made 
zero,  the  concentrated  transient  will  tend  to  have  a higher  maximum  crest 
than  maximum  trough.  If  the  constant  is  ± 0.5  the  maximum  crests  and 
troughs  are  the  same.  For  present  purposes  a value  of  -0.5  was  taken. 

All  that  remains  to  choose  in  the  evaluation  of  Eq.  (II)  is  the 
arbitrary  scaling  constant,  P,  in  the  specified  wavemaker  spectrum  modulus. 
This  constant  is  essentially  the  wavemaker  gain,  and  may  be  assigned  any 
number  from  zero  to  a value  at  which  the  waves  produced  will  be  too  steep 
and  start  to  break. 

Equation  (11)  was  evaluated  with  an  arbitrary  choice  of  P,  and  the 
parameters  selected  previously,  for  four  tank  positions  (at  the  wavemaker, 
and  at  125;  150  and  175  feet).  The  magnitude  of  the  resulting  wave  ampli- 
tudes was  compared  with  the  wave  length  corresponding  to  the  local  apparent 
periods  to  establish  a reasonable  maximum  value  of  P. 

Figure  3 shows  the  resulting  computed  transients  scaled  to  corres- 
pond to  the  estimated  maximum  amplitude.  As  required,  the  wave  pulse  at 
the  wavemaker  is  zero  for  time  less  than  zero.  The  time  duration  of  the 
pulse  decreases  with  position  in  the  tank  until  the  150  foot  position, 
whereafter  the  duration  increases.  The  transient  at  the  point  of  concen- 
tration ( 1 50 1 ) appears  small  in  the  figure  because  the  vertical  scale  is 
less  than  half  of  that  for  the  other  positions.  The  maximum  peak  to  peak 
amplitude  shown  for  the  150  foot  position  is  about  0.5  feet,  which  is 
probably  equivalent  to  a wave  height  to  wave  length  ratio  of  1/10.  At  the 
125  and  175  foot  positions  maximum  apparent  local  steepness  appears  to  be 
about  1/18,  and  at  the  wavemaker  maximum  apparent  local  steepness  is  less 
than  1/20.  For  these  reasons  the  wave  pulses  shown  in  Figure  3 were  taken 
as  representing  a practical  possibility  for  a relatively  severe  wave  pulse. 
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Having  established  that  the  computed  results  were  reasonable  for 
specific  points  in  the  tank,  some  examples  of  encountered  transients  could 
be  computed  from  Eq.  (13).  All  parameters  except  the  model  arrival  time 
delay,  6,  had  been  established,  and  Eq.  (13)  was  evaluated  for  a model 
speed  corresponding  to  a Froude  Number  of  0.15.  Three  cases  were  considered 
according  to  whether  the  model  arrives  at  the  150  foot  position  exactly  at 
time  t when  the  wave  pulse  is  most  concentrated  (6=0),  10  seconds  earlier, 
or  10  seconds  later. 

The  results  of  the  evaluation  are  shown  in  Figure  4 for  the  three 
cases.  As  in  the  previous  figure  the  vertical  scale  chosen  for  the  case 
of  model  arrival  during  maximum  concentration  is  smaller  than  that  for  the 
other  cases.  The  duration  of  the  encountered  transients  is  about  10  seconds 
in  all  cases.  This  duration  corresponds  to  roughly  4 model  lengths  travel 
up  the  tank.  For  the  speed  involved,  timing  of  model  arrival  at  a specific 
point  in  the  tank  within  plus  or  minus  a few  seconds  is  quite  feasible. 

The  complex  spectrum  was  computed  for  each  of  the  wave  pulses  shown 
in  Figure  4.  The  modulus  of  the  spectrum  of  the  three  transients  was  the 
same,  and  as  expected,  was  constant  within  the  anticipated  encounter  fre- 
quency range. 
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INPUT-OUTPUT  THEORY  FOR  RESISTANCE 

As  noted  in  the  introduction,  the  basic  assumption  in  the  present 
work  is  that  ship  model  resistance  follows  the  functional  polynomial 

JU 

input-output  model  treated  in  Refs.  1,2, 6, 9 , and  elsewhere.  This  approach 
involves  the  assumption  that  the  model  proceeds  at  constant  speed  in  a wave 
system,  h(t),  which  is  defined  at  a point  stationary  with  respect  to  the 
mean  position  of  the  model;  that  is,  the  input-output  model  applies  to  the 
encounter  domain. 

The  general  time  domain  representation  for  the  resistance,  r(t),  is 
wr i tten  as  fol lows : 

r(0  = 90  + «f9i  (t1)h(t-t1)dt1 

+ J,j92(ti>t2)h(t-ti)h(t“t2)dtldt2  (18) 

where: 

r(t)  is  instantaneous  resistance 

h(t)  represents  a zero  mean  wave  elevation 

t is  time 

and  gQ,  91(t1),  and  9a(tl,t2)  are  zeroth,  first  and  second  degree  kernels. 

It  is  presumed  that  the  hydrodynamic  properties  of  the  ship  model  are 
contained  wholly  in  the  kernels,  and  (without  loss  in  generality)  that 
the  second  degree  kernel  is  symmetric  in  its  arguments;  that  is, 

92(ta>ti)  = gaCt^tg) 

The  zeroth  degree  kernel,  gQ,  is  the  value  of  r(t),  Eq.  (18),  when 
the  wave  elevation  is  zero  for  all  time.  Since  it  may  reasonably  be  sup- 
posed on  physical  grounds  that  the  wave  elevation  will  be  zero  mean,  the 
kernel,  gQ,  is  identified  for  present  purposes  as  the  calm  water  resistance. 


9.  Dalzell,  J.F.,  "The  Applicability  of  the  Functional  Polynomial  Input- 
Output  Model  to  Ship  Resistance  in  Waves,"  SIT-DL-75-1 794,  Davidson 
Laboratory,  Stevens  Institute  of  Technology,  January  1975. 
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The  second  term  in  Eq.  (18)  is  a linear  convolution  integral  and 
is  interpreted  as  the  oscillatory  su^ge  exciting  force  so  that  the  first 
degree  kernel  is  the  impulse  response  function  for  surge  force. 

The  third  term  in  Eq.  (18)  is  a double  convolution  integral  which 
gives  rise  to  the  essentially  non-linear  resistance  added  by  waves.  The 
second  degree  kernel  may  be  called  a quadratic  impulse  response. 

Both  first  and  second  degree  kernels  are  assumed  to  be  absolutely 
integrable  and  thus  may  be  transformed  in  the  Fourier  sense.  The  trans- 
form pairs  relating  the  linear  and  quadratic  impulse  responses  to  corres- 
ponding linear  and  quadratic  frequency  response  functions  may  be  defined 
as  follows: 

9i  = hf  Je+IUJT  Gx  (u,)d(U 

Gj  (iu)  = Je‘lU>T  gx  (t)  dr  (19) 

MvV  = TT2  NExp  [+iVi  + iu^ 

(2rr)  u J 

= JJExP  r-iuJxTi  - (2°) 

in  which  a)  represents  circular  encounter  frequency. 

The  linear  frequency  response  function,  G1(tu),  defined  by  Eq.  (19) 
is  absolutely  conventional. 

The  quadratic  frequency  response  function,  Gs  (0^,11)  ) is  defined  in 
a bi  -frequency  plane.  Because  the  kernel  92(T.l>Ta)  is  assumed  to  be 
symmetrical  in  its  arguments,  and  is  real: 

GjKjWj)  = Gs  (21) 

G^,^)  = G2(-u,i>"u,3) 

= G2  (-U>a , -u^ ) (22) 

(The  asterisk  denotes  the  complex  conjugate.)  These  rel at ionsh i ps  simpli- 
fy the  quadratic  frequency  response  function  to  the  extent  that  certain 
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symmetries  result  and  that  as  a consequence  the  function  needs  only  to 
be  considered  in  a quadrant  of  the  bi-frequency  (<o  , u)a)  plane.  Equa- 
tion (21)  results  in  a line  of  symmetry  along  the  line  0^=0^.  Equation 
(22)  results  in  a line  of  symmetry  of  the  real  part,  and  ant i -symmetry 
of  the  imaginary  part  of  G„(oji,od2)  defined  by  =-uu1 . (It  may  be  noted 
that  along  this  line  the  imaginary  part  of  the  function  is  zero.)  These 
two  lines  and  the  (0^,(0^  axes  divide  the  bi-frequency  plane  into  octants, 
of  which  the  two  on  either  side  of  the  positive  axis  may  be  arbitrarily 
chosen  for  reference.  The  assumptions  of  symmetry  of  the  second  degree 
kernel  results,  with  Eq.  (20),  in  a complete  definition  of  G (uu,  ,u)  ) if 
the  functions  are  defined  in  any  pair  of  octants  including  a semi-axis 
of  either  frequency.  Thus  without  loss  in  generality,  interpretation  of 
the  quadratic  frequency  response  needs  only  to  involve  the  octants  on 
either  side  of  the  positive  axis.  In  these  octants  t»  is  positive  and 

! I * I I * 

Because  estimates  of  the  quadratic  frequency  response  function  are 
the  objective  of  the  present  work,  a summary  (after  Ref.l)  of  the  interp- 
retation of  the  function  is  in  order.  The  approach  to  the  meaning  of  the 
function  is  grossly  the  same  as  for  the  linear  case.  If  in  the  linear 
case  the  system  is  considered  to  be  excited  by 

h(t)  = a coscot 

the  output  may  then  be  written: 

Re  jaGx  (co)  Exp(  i u»t)  | 

and  G1(co)  is  interpreted  in  terms  of  normalized  amplitude  and  phase  of 
response. 

To  interpret  the  quadratic  frequency  response,  dua  1 harmonic 
excitation  is  necessary.  Accordingly,  it  may  be  assumed  that: 

h(t)  = ajCosu^t  + a costo  t (23) 

In  accordance  with  the  previous  discussion  of  symmetry,  both  frequencies 
(u^  ,<ua)  are  considered  positive  and  | J £ ] ui,  | . The  basic  model,  Eq.  (18) 
is  good  for  any  zero-mean  excitation.  Accordingly,  Eq.  (23)  may  be 
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substituted  directly  in  Eq.  ( 1 8) . After  some  algebra  the  final  result 
for  the  response  to  dual  harmonic  excitation  may  be  written  as  follows: 

= 90 

+ Re|ai  Gx  (u^  ) Exp(  iu^  0+a^  (u?2 ) Exp(  i t)  } 

+ (w1>-u)1)+ajG2(u)2,-u)s)  | 

+ | Reja^Gs  ,<1^)  Exp(i2u)x  t)  j 
+ | Reja^Gs(a)2,u)a)  Exp(i2iu2t)  | 

+ Reja1a2Gs(u)1,uia)Exp[i(tt)1+u)2)t]} 

+ Re|a1a2Ga((Ux,-(u?)Exp[i(aJ1-uj2)t]|  (24) 

This  result  shows  that  the  response  of  the  quadratic  system, 

Eq.  ( 1 8) , to  dual  excitation  contains,  in  general,  a shift  in  the  mean 
and  components  of  six  different  frequencies  [u^  , ujo  , 2^,  2uia , (tu^+u^), 
and  (ujj-U)  )].  The  second  and  third  terms  of  the  result  are  the  super- 
position of  the  linear  responses  at  the  excitation  frequencies.  The 
fourth  and  fifth  terms  of  Eq.  (24)  represent  a shift  in  the  mean.  These 
terms  allow  the  identification  of  the  mean  added  resistance  operator  as 
the  value  of  Ga  ( tUj , u)  ) along  the  line  0^=-^  (or  G2  (u^  , -u>  ) ) . The  sixth 
and  seventh  terms  are  the  second  harmonic  components  (2u^,2u}s).  Similarly, 
these  terms  allow  the  identification  of  second  harmonic  response  with  the 
values  of  G2(u>T,u)2)  along  the  line  ui  -Wj  (or  Ga  (u)^  ,u>l ) ) . 

The  eighth  and  ninth  terms  of  Eq.  (24)  pertain  to  the  bi -frequency 
plane  in  general.  The  eighth  term  is  the  response  at  frequency  (tUj+tu  ) ; 
that  is,  G (u\,ut,)  expresses  the  normalized  response  in  the  sum  frequency 
due  to  non-linear  interact  ions . Similarly,  the  ninth  term  involves 
response  at  frequency  (u^  -u)2)  ; that  is,  Gp  (u^  ,-ittj)  is  the  normalized 
response  in  the  difference  frequency. 

A final  part  of  the  general  theory  which  is  useful  in  both  analysis 
and  interpretation  is  the  effect  of  cascaded  linear  systems.  If  it  is 
supposed  that  the  input  to  the  system,  h(t)  of  Eq.  (18)  is  the  result  of 
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a linear  operation  upon  a variable  s(t),  and  the  frequency  response 
function  of  the  linear  operation  is  denoted  by  L (u>) , then  the  linear 


and  quadratic  frequency  response  functions  (K^  (uo)  and  (u^ 

ing  s(t)  and  r(t)  may  be  written: 

,u>a))  connect- 

(oj)  = G:  (uj)  (cu) 

(25) 

= G2('JJ1^a)L1  (<VLi  (“a) 

(26) 

If  it  is  additionally  supposed  that  the  output  r(t)  of  the  system  of 
Eq.  (18)  acts  upon  a linear  system  with  frequency  response  function  La  (uo) 
to  produce  an  output  v(t),  then  the  linear  and  quadratic  frequency  response 
functions  (K^V(uj)  and  K^V(u^  ) connecting  s(t)  and  v(t)  may  be  written: 

K^V(oj)  = (uj)  1^  (uj)  L2  (u>) 

(27) 

k|V(uu,,ui,)  = G2(ou1,uj2)L1  (uj1)L1  (ujs)L2(ui1+ujs) 

(28) 

The  relation  for  the  linear  part,  Eq.  (27),  is  just  the  same  as 
for  completely  linear  systems. 

Equation  (28)  shows  that  filtering  of  input  and  output  have  differ- 
ent effects  on  the  apparent  frequency  response  function.  For  example, 

> 

if  Lx  ( uj)  is  an  ideal  low-pass  filter  with  cutoff  frequency  a,  the  effective 
quadratic  frequency  response  function,  K^V(uj1  ,u;2) , will  be  zero  outside  the 
region  where  I j and  | u)2  j are  less  than  <?.  Considering  a practical 
example  of  output  filtering,  if  the  output  is  observed  through  a low-pass 
filter,  L2  (ui)  , with  a very  low  frequency  cutoff,  the  effective  quadratic 
frequency  response  will  tend  to  be  zero  everywhere  except  near  the  line 

which  is  the  location  of  the  mean  added  resistance  operator  in  the 
bi-frequency  plane. 
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SIMULATION  OF  RESISTANCE  TRANSIENTS 


Genera  I 

Following  the  conditions  involved  in  the  simulation  of  encountered 
wave  pulses.  Figure  4,  it  was  assumed  that  the  simulated  experiment 
involved  a five  foot  ship  model  proceeding  at  a Froude  Number  of  0.15  into 
a head  wave  pulse.  In  order  to  take  advantage  of  quantitative  resistance 
response  data  developed  in  Refs,  2,  6 and  9,  the  model  was  assumed  to  be 
the  Series  60  0.60  block  parent,  restrained  in  surge.  Because  it  is 

customary  in  a narrow  towing  tank  to  locate  a wave  probe  ahead  of  the  model 
to  minimize  distortion  of  results  by  model  generated  waves,  the  encountered 
wave  pulses,  Figure  4,  were  assumed  to  be  the  indication  of  a wave  probe 
one  model  length  forward  of  model  LCG. 

The  simulated  wave  pulses  are  digital  time  series;  that  is, 
represent  a sampling  of  the  wave  pulse  at  a uniform  time  interval  (At). 
Accordingly  the  analytical  time  domain  model,  Eq.  (18),  has  to  be  re-cast 
into  a summation  form  for  practical  computation  purposes.  This  was  done 
in  the  same  way  as  had  been  done  in  Ref.  9.  The  digital  model  of  Eq.  (18) 
becomes : 

r(n)  = gQ 

m2 

+ EL.  h(n-j) 

j =-m  1 ■* 

p2  P2 

+ E E Q h(n-j ) h(n-k)  (29) 

j =—  p 1 k=-p  1 -1 

where  r(n)  = computed  resistance  time  series 
h(n)  = simulated  wave  pulse  time  series 

and  L.  and  Q...  are  weighting  coefficients  in  which  the  differentials  are 
J J K 

absorbed . 
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Equation  (29)  requires  discrete  approximations  to  the  linear  and 
quadratic  impulse  responses  of  the  form: 


gt (t)  = L.  6 (r-jAt) 

(30) 

9s(ti,ts)  = Qjk  6(Tx-jAt)6(T2-kAt) 
(where  6(t)  is  the  Dirac  delta  function) 

(31) 

The  Fourier  transforms  of  Eqs.  (30)  and  (31)  in  accordance  with  Eqs.  (19) 
and  (20)  yield  continuous  (aliased)  frequency  response  functions  as 
fol lows : 

m2 

^ (uj)  = E L . [cos  (j uAt)  - i sin(juAt)]  (32) 

j =-ml  J 

p2  P2 

G_(ul,<b)  = S E Q.,  [cosQu^At+kiu  At)-i  s i n(ju) At+kui At)  ] (33) 

* ^ ® • ii  i I K 2 * Z 

J =-p 1 k=-p I J 


For  practical  computing  purposes  the  process  of  transforming  given 
observed  or  analytical  frequency  response  functions  amounts  to  achieving 
(with  a finite  number  of  coefficients)  a reasonable  fit  of  the  re-trans- 
formed impulse  responses  (Eqs.  32,33)  to  the  observed  response  functions. 
The  general  approach  followed  is  to  integrate  the  first  of  Eqs.  (19)  and 
(20)  trapezoidal ly  with  observed  estimates  of  the  appropriate  frequency 
response  function.  This  is  carried  out  for  a sufficient  range  of  time 
variable(s)  so  that  decisions  regarding  the  truncation  limits  (ml,  m2, 
and  pi,  p2  in  Eqs.  29  through  33)  can  be  made.  Evaluation  of  the  computed 
impulse  response  function  is  made  at  integer  values  of  time  step  At  and 
these  values  are  multiplied  by  At  or  At  as  appropriate  to  result  in  esti- 
mates of  the  coefficients  L.  and  Q..  . The  next  step  is  to  adjust  the 
coefficients  so  that  the  discrete  kernel  reflects  the  correct  behavior 
of  the  frequency  response  functions  at  zero  frequency.  In  the  present 
case  it  is  assumed  that  there  is  no  resistance  added  by  an  infinitely 

long  wave.  Thus  the  coefficients  L.  and  Q.  should  sum  to  zero.  The 

J J k 

final  computational  step  is  to  insert  the  coefficients  into  Eqs.  (32)  or 
(33),  compute  the  re-transformed  impulse  response  function  and  compare 
the  results  with  the  original  frequency  response  estimates.  Iteration 
of  the  process  is  sometimes  necessary. 
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Selection  of  Data  for  Use  in  Eg.  (59) 

To  take  advantage  of  the  results  of  Refs.  2,  6 and  9,  it  was 
convenient  to  retain  the  non-d imens ional izat ions  of  those  references. 

The  first  convention  is  that  all  lengths  are  divided  by  model  length. 

For  present  purposes  this  meant  only  that  the  simulated  wave  pulse  eleva- 
tions were  to  be  divided  by  five  feet.  The  second  convention  was  that 
model  resistance  was  to  be  divided  by  model  displacement.  Finally, 
a non-dimensional  encounter  frequency  was  defined  as: 

o = ou/  0)^  L 

where  tu  = encounter  frequency,  rad/sec 

= frequency  of  a wave  of  model  length,  L 

= /2^7l 


Because  this  non-d imens i ona 1 i zat ion  can  be  considered  as  merely  a 
change  in  the  time  scale,  the  a notation  of  Refs.  2,  6 and  9 may  be  used 
interchangeably  with  the  present  ui  notation.  For  present  purposes  the 
time  scale  will  be  multiplied  by  U)^  and  corresponding  frequencies  com- 
puted in  the  conventional  way.  In  particular,  the  real  time  sampling 
interval  (At)  of  the  wave  pulses  of  Figure  k was  selected  to  be 
0.066  seconds.  The  equivalent  At  in  the  "non-dimensional"  time  scale  is 
0.4197>  and  this  value  was  used  through  the  developments  of  the  weighting 
coefficients  in  Eqs.  (29)  through  (33)  to  be  described. 


The  first  term  in  Eq.  (29)  is  the  calm  water  resistance,  g^.  The 
experimental  value  of  raw  calm  water  resistance  cited  in  Ref.  6 for  the 
five  foot  series  60  model  at  a Froude  Number  of  0.15  was  0.092  lb. 
Dividing  this  by  model  displacement  yields: 

g = 0.0028 
3o 

The  second  term  in  Eq.  (29)  involves  the  linear  weighting  coef- 
ficients L..  Several  sets  of  estimates  from  experiment  had  been  made  in 
Ref.  6 for  the  modulus  of  the  surge  exciting  force  frequency  response 
function,  Gx  (u>)  or  GT  (a) . These  were  averaged  to  form  a starting  point 
for  the  derivation.  The  phase  data  was  not  presented  in  Ref.  6,  but 
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had  been  developed.  These  phase  estimates  had  been  derived  with  respect 
to  wave  elevation  forward  of  LCG.  This  wave  reference  had  been  shifted 
in  time  to  improve  the  cross-spectral  estimation  procedure,  and  it  was 
a straightforward  operation  to  correct  the  estimates  to  a wave  phase 
reference  one  model  length  ahead  of  model  LCG. 

The  operations  described  in  conjunction  with  Eqs . (30)  through 
(33)  were  carried  out  upon  the  experimental  surge  exciting  force  data 
and  the  results  are  summarized  in  Figures  5 and  6.  Figure  5 indicates 
the  experimental  estimates  of  the  surge  exciting  force  frequency  response 
function  as  plus  signs,  and  the  results  of  re-transforming  the  final 
weighting  coefficients,  L^,  as  dashed  lines.  As  may  be  noted,  the 
experimental  data  which  was  thought  valid  in  Ref.  6 does  not  cover  the 
entire  frequency  range.  (Extrapolations  had  been  made  to  zero  and  high 
frequency  before  starting  the  transformation  procedure.)  However  the 
encounter  frequency  range  of  the  example  wave  pulses  does  not  significantly 
exceed  the  frequency  range  where  data  is  shown,  and  it  would  thus  not  be 
expected  that  the  results  of  the  simulation  would  be  seriously  influenced 
by  errors  in  the  extrapolations.  Overall,  the  correspondence  between  the 
transformed  discrete  linear  kernal  and  the  experimental  data  appears  ade- 
quate to  insure  realism  in  the  simulation. 

The  weighting  coefficients,  L^,  in  Figure  6 are  of  most  importance 
for  positive  k;  that  is,  for  the  "past"  of  the  input  wave  elevation. 
Accordingly  the  kernel  relating  surge  exciting  force  to  wave  elevation 
a model  length  ahead  of  LCG  is  nearly  physically  realizable. 

For  a starting  point  for  the  simulation  of  the  third  term  in 
Eq.  (29)  there  were  four  choices  available.  Reference  9 contains  three 
sets  of  experimental  estimates  of  the  quadratic  frequency  response  function 
for  the  experimental  case  of  interest,  and  Ref.  2 contains  one  set  of 
analytical  estimates.  It  appears  from  the  analyses  of  Ref.  2 that  in  the 
regions  of  the  bi-frequency  plane  not  resolved  in  Ref.  9 there  exist 
significant  values  of  the  quadratic  frequency  response  function.  Addi- 
tionally, some  of  the  experimental  estimates  which  were  questioned  in 
Ref.  9 appear  to  have  been  discredited  by  the  analyses  of  Ref.  2. 
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ENCOUNTER  FREQUENCY,  <J 


FIGURE  5 SURGE  EXCITING  FORCE  FREQUENCY  RESPONSE  FUNCTION 

FOR  SERIES  60,  0.60  BLOCK  MODEL  AT  FROUDE  NUMBER  0.15 
(PHASES  ARE  REFERENCED  TO  WAVE  ELEVATION  ONE 
MODEL  LENGTH  FORWARD  OF  LCG) 
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FIGURE  6 COEFFICIENTS,  lk,  OF  DISCRETE  LINEAR 

IMPULSE  RESPONSE  CORRESPONDING  TO 
FREQUENCY  RESPONSE  FUNCTION  OF  FIGURE  5 
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However,  in  the  region  of  the  bi-frequency  plane  of  most  practical 
interest  (that  is,  near  the  mean  added  resistance  operator)  the  experi- 
mental and  analytical  estimates  are  not  impossibly  different.  Accord- 
ingly, the  best  choice  for  present  purposes  appeared  to  be  the  analytical 
result  of  Ref.  2. 


The  analytical  quadratic  frequency  response  function  of  Ref.  2 
relates  the  quadratic  part  of  added  resistance  to  wave  elevation  at  LCG. 
Thus  the  first  necessary  operation  for  the  simulation  was  to  develop  the 
analytical  quadratic  frequency  response  function  relating  resistance  to 
wave  elevations  one  model  length  ahead  of  LCG.  This  was  accomplished  in 
accordance  with  the  input  filtering  relationship,  Eq.  (26).  Essentially, 
wave  elevation  ahead  of  model  and  at  LCG  are  related  in  the  wave  frequency 
domain  by  the  transfer  function  following  Eq.  (3).  This  transfer  function 
becomes  La  (cu)  after  transformation  into  the  encounter  frequency  domain. 


Having  thus  made  an  initial  estimate  of  the  quadratic  frequency 

response  function,  the  process  outlined  in  the  discussion  of  Eq.  (33)  was 

pursued.  The  final  results  for  the  Q.  weighting  coefficients  correspond- 

J * 

i ng  to  the  quadratic  impulse  response  are  shown  in  an  isometric  view  in 
Figure  7.  The  weighting  coefficient  matrix  is  71  x 71,  and  the  value  of 
the  largest  coefficient  is  0.69.  In  the  figure  straight  lines  are  drawn 
through  the  points  along  sections  parallel  to  j and  k axes  so  that  each 
intersection  corresponds  to  a value  of  (J.^.  Where  the  function  is  sig- 
nificant it  is  character ized  by  undulations  in  both  j and  k directions. 
Most  of  the  negative  values  of  the  function  are  masked  or  indistinct, 
but  are  of  roughly  the  same  magnitude  as  the  positive  peaks  shown. 

As  may  be  noted  at  the  edges  of  the  plot,  the  truncation  could  have  been 
at  slightly  higher  values  of  j and  k,  but,  as  will  be  shown,  the  repre- 
sentation seems  adequate.  In  this  plot  as  in  Figure  6 positive  values 
of  j and  k correspond  to  the  "past"  of  the  wave  elevation.  The  most 
significant  part  of  the  kernel  lies  in  the  region  where  both  j and  k 
are  positive.  There  are  weak  interactions  between  positive  k and  nega- 
tive j (and  positive  j and  negative  k by  symmetry),  but  very  little 
interaction  response  when  both  j and  k are  negative. 
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Agreement  was  quite  good  between  the  quadratic  frequency  response 
function  corrected  to  wave  forward  of  LCG  and  the  re-transformation  of 
the  truncated,  discrete  function  shown  in  Figure  7.  The  re-transformed 
function  was  then  corrected  to  a wave  phase  reference  at  model  LCG  and  is 
shown  in  comparison  with  the  original  analytical  result  of  Ref.  2 in 
Figures  8 and  9.  Figure  8 indicates  the  real  part  and  Figure  9 the  imag- 
inary part  of  the  function.  These  two  figures  are  essentially  plotted 
tables.  The  magnitude  of  the  function  involved  is  plotted  in  the  non- 
dimensional  bi-frequency  plane  (0^,11^).  A second  axis  system  is  shown  as 
well.  This  axis  system  is  the  "sum"  and  " difference " system  defined  by 
the  transformations: 

Ol  = “i-Ufe 

Qs  = ujj+ut,  (34) 

The  and  axes  lie  on  the  lines  of  symmetry. 

The  results  from  the  re-transformation  of  the  discrete  kernel, 

Figure  7,  are  shown  in  the  quadrant  bounded  by  the  +Q1  and  Qs  axes.  In 
order  to  show  the  analytical  results  from  Ref.  2 in  the  same  figures,  the 
octant  bounded  by  the  +iu2  and  +QS  axes  has  been  displaced  in  the  negative 
C3X  direction.  The  analytical  results  are  shown  in  this  octant  according 
to  the  basic  symmetry  properties  of  the  quadratic  frequency  response 
function.  If  analytical  results  and  the  simulation  were  in  exact  agree- 
ment, the  simulated  results  would  be  an  exact  reflection  about  the  axis 
of  the  analytical  results  in  the  uj  -Qa  octant.  A similar  procedure  was 
adopted  to  show  the  analytical  results  corresponding  to  the  simulated 
results  in  the  octant.  In  this  case  exact  agreement  would  require 

that  the  simulated  results  would  be  an  exact  reflection  about  the  C3  axis 
of  the  analytical  results  in  the  ^-(-u^)  octant. 

In  order  to  emphasize  the  important  parts  of  the  function,  esti- 
mates having  absolute  values  less  than  2.5  are  omitted  (blank).  The 
contours  shown  are  for  the  zero  level. 
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ANALYTICAL 
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On  the  whole  the  correspondence  between  simulation  and  analytical 
result  is  considered  quite  adequate  for  present  purposes.  The  re-trans- 
formation was  carried  out  at  half  the  frequency  intervals  shown  in  the 
figures.  The  intermediate  values  indicated  that  the  function  was  smooth. 

Computation  of  Sample  Transient  Responses 

Having  selected  coefficients  of  Eq.  (29)  which  were  thought  to 
fairly  realistically  represent  the  resistance  characteristics  of  the 
assumed  ship  model,  it  remained  only  to  do  the  convolution  arithmetic 
for  the  example  encountered  wave  pulses  of  Figure  4. 

There  is  a numerical  start-up  and  stopping  transient  in  both 
convolutions  in  Eq.  (29).  However  in  the  present  case  the  input  wave 
pulse  is  zero  at  beginning  and  end  so  that  the  resulting  outputs  do  not 
have  to  be  specially  treated  for  this  effect. 

The  basic  programming  used  was  that  developed  in  Ref.  9.  Before 
computing  the  response  to  the  wave  pulses  the  programming  was  run  with 
sinusoidal  wave  elevation  input,  and  analysis  of  results  in  accordance 
with  Eq.  (24)  confirmed  that  the  computing  system  was  correct. 

As  a possible  aid  in  interpretation,  the  components  of  Eq.  (29) 
were  separately  evaluated  for  the  wave  pulses  of  maximum  amplitude  (Fig. 4) 
and  these  results  were  stored  separately  for  later  use,  so  that  the 
results  of  the  computation  were  essentially  three  1024  point  time  series: 

1.  The  non-dimensional  input  wave  pulse  of  maximum 

ampl i tude 

2.  The  resulting  linear  component  of  Eq.  (29) 

3.  The  resulting  quadratic  component  of  Eq.  (29) 

To  produce  a simulated  resistance  transient  which  could  be  considered  as 

the  response  to  a wave  pulse  of  amplitude  (F)  times  the  maximum  shown  in 

Figure  4,  the  linear  component  time  series  is  multiplied  by  F,  the  quad- 

2 

ratic  component  time  series  is  multiplied  by  F , and  the  simulated 

transient  is  the  sum  of  these  modified  series  and  q . 

o 

I 
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The  basic  results  of  the  computation  are  shown  in  Figures  10 
through  13.  The  wave  pulses  of  Figure  4 were  numbered;  No.  1 corres- 
ponds to  model  arrival  at  mid-run  10  seconds  early,  No.  2 to  arrival 
on  time,  and  No.  3 to  model  arrival  late.  Figure  10  involves  wave  pulse 

No.  1,  Figure  11  pulse  No.  2 and  Figure  12  pulse  No.  3 

In  each  of  the  figures  the  non-d imens ional ized  wave  pulse  is  shown 
in  the  top  frame,  the  linear  and  quadratic  components  of  resistance  in  the 
next  two  frames,  and  finally  at  the  bottom  the  simulated  resistance  (sum 
of  the  calm  water  resistance  and  the  linear  and  quadratic  components). 

The  time  scale  shown  in  the  figures  is  essentially  the  point  number 
of  the  time  series.  Only  300  points  of  the  1024  available  were  plotted 
for  each  case.  The  portions  of  the  time  series  not  shown  are  constant 
and  equal  to  the  values  at  beginning  and  end  of  the  portions  plotted. 

The  vertical  scales  for  each  time  history  frame  were  chosen  to  best 
resolve  the  transient  being  plotted  and  are  not  the  same  from  frame  to 
frame. 

Figure  13  illustrates  the  simulated  results  when  the  wave  pulse 

is  assumed  to  be  less  steep  than  maximum.  In  this  case  the  simulated 

results  for  wave  pulse  No.  1 of  maximum  (full)  amplitude  are  compared 
with  the  simulations  for  the  cases  when  the  wave  pulse  is  assumed  to 
involve  wave  elevations  1/2  and  1/4  the  maximum  shown  in  Figure  4. 
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DISCUSSION  OF  RESULTS  OF  RESISTANCE  SIMULATION 


In  review  of  the  simulated  results  presented  in  previous  sections 
it  may  first  be  repeated  that  most  of  the  detail  was  undertaken  in  an 
attempt  to  provide  realistic  results  for  one  particular  choice  of  model 
and  experimental  parameters.  The  experimental  observations  are  concept- 
ually confined  to  the  wave  pulse  as  observed  by  a wave  probe  one  model 
length  ahead  of  the  model  (top  frames  of  Figures  10  through  13)  and  the 
resul ti ng  'res i stance  transient  (bottom  frame  of  Figures  10  through  12). 
These  simulated  results  are  defined  as  time  series  in  the  same  way  as 
actual  experimental  data  would  be,  after  digitization.  The  components  of 
resistance  transients  shown  in  Figures  10  through  12  are  not  observable. 
The  frequency  response  functions  shown  in  the  previous  section  are  also 
not  available  — these  (or  parts  thereof)  are  the  answers  sought  in  the 
analysis  of  the  transients. 

Within  the  particular  assumptions  of  the  hypothetical  experiment 
(O.L.  Tank  No.  3,  a five  foot  Series  60  model  at  Froude  Number  0.15)  the 
wave  pulses  developed  appear  to  be  realistic  possibilities,  and  the  total 
length  of  the  simulated  runs  (including  portions  not  plotted  in  Figures 
10  through  13)  is  feasible.  It  would  also  be  feasible  to  produce  wave 
pulses  of  longer  duration  by  altering  the  assumptions  of  a previous 
section  with  respect  to  the  position  in  the  tank  of  maximum  wave  concen- 
tration. 

Given  the  feasibility  of  obtaining  relatively  long  stretches  of 
essentially  calm  water  resistance  at  beginning  and  end  of  the  observed 
resistance  transient,  the  feasibility  of  the  first  and  most  obvious 
step  in  the  ans lysis  of  the  resistance  transient  is  confirmed.  This  is 
to  remove  the  effect  of  calm  water  resistance  (gQ)  by  fitting  a straight 
line  through  initial  and  final  portions  of  the  observed  transient,  and 
then  subtract  this  mean  line  from  the  observation.  The  portion  of 
observation  left  after  this  correction  should  be  the  sum  of  linear  and 
quadratic  components  due  to  wave  excitation.  In  practice,  since 
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electronic  zero  bias  and  possibly  slight  trends  must  be  dealt  with, 
this  type  of  procedure  would  probably  have  to  be  carried  out  for  the 
wave  pulse  as  well  so  as  to  produce  "zero"  wave  elevation  at  beginning 
and  end  of  the  observations. 

The  individual  components  of  resistance  shown  in  Figures  10 
through  12  are  about  what  would  be  expected  from  the  input-output  model. 

The  linear  component  of  resistance  resembles  a time  lagged  wave  pulse 
insofar  as  symmetry  about  zero  and  apparent  periods  are  concerned. 

The  quadratic  component  does  not  much  resemble  the  wave  pulse  except 
in  duration.  This  component  contains  significant  response  at  about  half 
the  apparent  periods  in  the  wave  pulse  and  a visible  relatively  long 
period  excursion  of  resistance  so  that  the  quadratic  component  of  resist- 
ance is  not  symmetric  about  the  mean. 

The  simulated  total  resistance  transients  in  Figures  10  through  12 
are  visibly  influenced  by  the  quadratic  component.  These  results  are  for 
a wave  pulse  of  probable  maximum  intensity.  It  may  be  noted  from  Figure  13 
that  the  appearance  of  the  resistance  transient  is  significantly  influ- 
enced by  the  magnitude  of  the  wave  pulse.  A comparison  of  the  linear 
component  of  Figure  10  with  the  resistance  for  a wave  pulse  of  1/4 
maximum  amplitude  in  Figure  13  discloses  that  if  the  gain  of  the  hypothet- 
ical wavemaker  is  reduced  so  that  the  wave  steepnesses  are  1/4  as  great 
as  the  probable  maximum  obtainable,  the  obvious  influence  of  the  quadratic 
non-linearity  on  the  total  observed  resistance  transient  all  but  dis- 
appears . 

The  foregoing  is  of  course  implicit  in  the  input-output  model 
assumed,  and  was  observable  in  a sense  in  the  irregular  wave  experiments 
described  in  References  6 and  9.  In  the  examples  herein  the  quadratic 
component  of  resistance  is  of  the  same  magnitude  as  the  linear  part  only 
for  the  wave  pulses  of  an  amplitude  near  the  limit  for  propagation  with- 
out substantial  change  in  form.  In  practice,  with  recording  equipment  of 
more  or  less  fixed  dynamic  range,  it  would  be  somewhat  natural  to  opt  for 
a relatively  severe  wave  pulse  when  best  resolution  of  the  quadratic  com- 
ponent is  desired.  On  the  other  hand,  the  qualitative  indications  of 
Figure  13  suggest  that  isolation  of  the  linear  component  might  be  achieved 
to  acceptable  accuracy  if  wave  pulses  of  the  mildest  amplitude  possible 
are  employed. 
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INITIAL  ANALYSIS  OF  QUADRATIC  TRANSIENT  RESPONSE 


Existing  Theoretical  Methods 

Among  the  most  fundamental  and  useful  engineering  treatments  of 
the  functional  input-output  model  (Barrett  ^*  , Bedrosian  and  Rice'*  , 

George  ) only  George  treats  the  possibility  of  recovering  the  quad- 
ratic impulse  or  frequency  response  from  transient  experiments.  The 
approach  of  the  latter  author  involves  successive  isolation  of  the  time 
domain  kernels  in  a series  of  experiments.  The  first  "experiment"  is 
essentially  the  same  as  that  implied  in  previous  discussion  and  consists 
of  observing  that  the  system  output  for  zero  input  is  the  zeroth  degree 
kernel,  gQ.  The  second  series  of  "experiments"  involves  excitation  of 
the  system  by  ideal  step  functions  of  varying  amplitudes.  The  resulting 
sets  of  output  transients  are  expanded  in  a series  in  the  amplitude  of 
the  step  function  for  each  of  a number  of  successive  values  of  time. 

The  coefficients  of  the  first  power  of  the  step  amplitude  form  a function 
of  time  and  the  linear  impulse  response  is  its  derivative  with  respect  to 
time.  The  third  series  of  "experiments"  involves  excitation  of  the  system 
with  two  step  functions  which  are  separated  by  a known  time  delay,  t. 

This  is  repeated  for  a series  of  step  amplitudes  as  before,  and  a series 
expansion  in  step  amplitude  is  made.  After  utilizing  some  results  of  the 
second  experiment,  the  coefficients  of  the  second  power  of  the  step  ampli- 
tude form  a function  of  time,  t,  and  time  delay,  t.  The  whole  procedure 
must  be  repeated  for  different  values  of  t so  that  finally  a two-dimensional 
function  of  time  (say  Tj  and  t ) is  generated.  The  quadratic  impulse 
response  is  effectively  the  partial  derivative  with  respect  to  and  t . 
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The  above  procedure  depends  critically  upon  the  properties  of 
ideal  step  functions,  and,  accordingly,  this  procedure  is  unworkable 
for  the  present  problem.  It  is  not  possible  to  generate  wave  steps, 
ideal  or  otherwise.  In  fact,  judging  by  the  numerical  results  shown  in 
Figure  4,  it  is  doubtful  that  it  would  be  feasible  to  so  control  the 
various  experimental  parameters  that  wave  pulses  of  exactly  the  same 
shape  but  different  amplitudes  could  be  encountered  by  the  model.  Even 
if  wave  steps  were  physically  possible,  the  multiple  run  requirement  would 
make  the  approach  unattractive  since  the  only  point  in  developing  a 
transient  test  technique  is  to  make  one  test  run  do  the  work  of  many. 


Correlation  and  Fourier  Integral  Analysis 

First  intuition  suggested  that  analysis  of  the  non-linear  trans- 
ient might  be  approached  in  grossly  the  same  manner  as  the  previous 
analyses  of  the  random  process  with  quadratic  non-linearity'.  Clearly, 
most  of  the  techniques  available  for  analysis  of  linear  random  systems 
have  their  counterpart  in  the  analysis  of  linear  transients.  The  most 
important  difference  is  that  the  operation  of  taking  the  statistical 
expectation  in  the  random  analysis  is  replaced  in  the  transient  analysis 
by  integration  over  all  time.  The  expectation  or  mean  value  of  the  type 
of  wave  induced  transients  shown  herein  clearly  tends  to  zero  as  the 
limits  of  integration  are  increased  without  bound. 


In  the  analyses  to  follow  the  n-d imens i onal  form  of  Parseval's 
Theorem  (Barrett'^)  was  used  extensively,  and,  for  reference,  it  is 
reproduced  as  follows: 


If’ 


’Jfi  (ti>  V 


tn>dti 


dt„ 


• dt 


(2nj 


n !!•  • • XK ’ “s  • • ‘X  Fs  H '“fe  ' " ‘X  du)i  du)2  • * ,du;n 


(35) 


where  the  (*)  denotes  complex  conjugate  and  f ^ ( t4  • * • ) and  F.(t1***)  are 
Fourier  Transform  pairs  defined: 
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n 

F 5 (u>,,U)s  •••(*)„)  = Jj-'vff  i(ti>ta*"tn)ExP("'  Z u>rt  )dt.dt  ...dt 

J J * r=l 

I n 

" 7~7Z  i’J’”-IFjK^3-**“,n)Exp(i  z1a,rtr)d“l",d“n  (36) 

J (2tt)  j r- 1 


The  simplest  operation  on  the  resistance  transient  is  integration 
over  time.  Using  the  basic  model,  Eq.  (18),  and  transposing  the  constant, 

^o* 

J(  r ( t) -gQ) dt  = JJgj (t1)h(t-t1)dtdtl 

+ Jj9a(t1,ta)Jh(t-t1)h(t-ta)dtdtldta  (37) 

Noting  that  the  first  term  on  the  right  hand  side  is  zero  for  a zero- 
mean  wave  pulse,  and  applying  the  Parseval  Theorem: 

j(r(t)-gQ)dt  = ^ jGs  (uj, -id)  | H(uj)  | S dtu  (38) 

where  H ( uu)  is  the  Fourier  Transform  of  the  wave  pulse.  The  integral  of 
the  wave  induced  part  of  r(t)  is  non-zero  in  general  and  in  form  is  the 
same  as  the  expression  for  the  expected  value  of  resistance  in  random 
seas',  the  squared  absolute  value  of  the  complex  wave  pulse  spectrum 
being  analogous  to  the  scalar  spectrum  of  the  random  waves.  The  operation 
yields  no  approach  to  the  identification  of  the  mean  added  resistance 
operator,  G2  (u), -id). 

The  next  integral  operation  which  is  common  in  analysis  of  the 
linear  case  is  a lagged  product.  Forming  the  lagged  product  of  resist- 
ance transient  and  wave  pulse,  and  substituting  Eq.  ( 1 8) : 

mi = Jr(t)h(t-x)dt 

= 9QJh(t-T)dt 

+ JJSi  (t1)h(t-tl)h(t-T)dtldt 

+ JJ,j98(t1,t2)h(t-t1)h(t-t8)h(t-T)dtldtadt  (39) 

The  first  observation  about  Eq.  (39)  is  that  the  integral  multiplied  by 
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gQ  is  zero  within  the  zero  mean  wave  pulse  assumption.  Utilizing  the 
Parseval  Theorem  to  modify  the  succeeding  integrals  there  results: 

1 r 

(uu)  jm1  (*)  Exp[-iuoi]dT 

= (up)  | H(oj)  1 ? 

* h "*<■)  (40) 

where  inside  the  second  integral  the  frequency  parameters  uj  and  (3  corres- 
pond to  the  sum  and  difference  frequencies,  Qs  and  which  are  defined 
by  Eqs.  (34).  If  the  system  had  been  entirely  linear  the  integral  would 
be  zero  and  the  result  would  be  exactly  as  expected  for  linear  systems. 

In  the  case  of  random  wave  excitation  the  analogous  integral  dropped  out 
because  the  excitation  was  assumed  to  be  zero-mean  Gaussian,  and  the 
expected  value  of  triple  products  of  Gaussian  variables  is  zero.  Thus 
while  simple  cross-correlation  analysis  is  fruitful  in  identification  of 
the  linear  term  from  random  data,  it  does  not  appear  useful  in  the  case  of 
transient  data  since  only  the  elimination  of  the  constant  term  is  effected. 

Proceeding  with  the  analogy  with  the  analysis  of  random  data, 
a double  lagged  product  may  be  defined: 

m2(VTs)  = Jh(t+T1)h(t-T1)[r(t-T2)-g0]dt 

= (t1)h(t-T2-t1)h(t+T1)h(t-Tl)dt1dt 

+ JIJ92(t1,t2)h(t-i:s-tl)h(t-Ta-ts)h(t+-tI)h(t-Tl)dt1dtrdt 

(41) 

This  product  is  essentially  the  same  as  the  product  used  in  the  deriva- 
tion of  the  cross-bi -spectrum.  The  difference  is  in  the  integration  with 
time  instead  of  the  expected  value  operator. 

Applying  the  Parseval  Theorem  to  the  integrals  of  Eq.  (41)  there 
resul ts : 
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MV")  = Jma<Ti*Tai  Exp[-iu)Ta]dTa 


where 


= H2(u),t1)jg“(uj)h"((u) 


H2(uj,t)  = Jh(t--r)  h(t+T)  Exp[- i tu)]dt 


(42) 

(43) 


and  a)  and  £}  of  the  second  integral  are  the  same  as  those  in  Eq.  (40). 

The  terms  within  curly  brackets  differ  from  Eq.  (40)  virtually  only  by 
a factor  of  H(uj)  . The  quadratic  frequency  response  function  is  imbedded 
in  an  integral.  No  separation  of  linear  and  quadratic  frequency  response 
functions  has  been  effected  for  exactly  the  same  reasons  as  in  the  previ- 
ous case  for  the  single  lagged  product.  The  formation  of  the  lagged 
product  as  a function  of  is  seen  to  be  redundant.  If  t1  is  considered 
constant  and  equal  to  zero  in  Eq.  (4 1 ) , the  function  of  which  results 
corresponds  to  the  lagged  product  of  resistance  and  the  squared  wave 
pulse.  The  result  of  this  variation  is  exactly  the  same  in  form  as 

Eq.  (42). 

A final  integral  analysis  of  interest  is  the  Fourier  Transform  of 
the  resistance  transient,  R(u))  . The  resistance  transient  and  its  Fourier 
transform  are  a transform  pair  as  follows: 


8(uj)  = Jr(t)  Exp[-io)t]dt 

(44) 

r(t)  = Jr(uj)  Exp[  iu)t]du> 

(45) 

Manipulating  the  integrals  in  the  basic  input-output  model,  Eq.  (18),  by 
means  of  the  Parseval  Theorem, the  Fourier  transform  of  the  resistance 
transient  becomes: 

R(<d)  = 2rrgQ  6 (a))  + G1(uu)H(uj) 

(-6) 

where  6 (uu)  is  the  Dirac  delta  function.  This  is  the  expression  analogous 
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to  the  scalar  spectrum  of  resistance  in  the  random  excitation  case  , and 
the  form  is  very  similar.  Since  H(0)  must  be  zero  for  a zero  mean  wave 
pulse,  R(0)  may  be  seen  to  be  the  integral  of  the  wave  induced  resistance, 
Eq . (38),  plus  a delta  function  arising  from  the  constant  calm  water 
res i stance. 

Reviewing  the  results  in  Eqs . (40) , (42)  and  (46)  it  is  clear  that 
the  correlation  and  Fourier  Integral  approaches  do  not  immediately  suggest 
any  general  methods  for  data  reduction.  No  separation  of  linear  and  quad- 
ratic terms  is  evident.  The  quadratic  frequency  response  function  appears 
in  an  integral,  and  this  integral  is  practically  the  same  regardless  of 
the  initial  approach.  Simple  Fourier  Transformation  of  the  resistance 
transient  and  wave  pulse  results  in  information  of  exactly  the  same  nature 
as  that  found  by  correlation  techniques. 

Some  explanation  of  the  result  is  afforded  by  supposing  that  the 
quadratic  impulse  response  is  known,  as  well  as  an  input  transient  h(t) . 
Under  these  circumstances  a function  of  two  times,  t and  t,  may  be 
developed : 

p(t>T)  = Jl92(t1»t2)h(t-t1)h('r-t2)dt1dt2  (47) 

Applying  the  Parseval  Theorem  to  Eq.  (47): 

P(oUi,u)s)  = jjp(t,x)  Exp[-ia)lt-i^T]dtdT 

= G2(u)1,u)s)H(u^)H(u)2)  (48) 

Given  the  function  p(t,x)  and  the  input  transient  h(t)  Eq . (48)  implies 
a simple  identification  method  for  G.,  (u^  ,u^)  . However,  it  may  be  seen 
by  comparing  Eq.  (47)  with  the  quadratic  term  of  Eq . ( 1 8)  that  the 
quadratic  component  of  resistance  is  the  value  of  the  function  p(t,x) 
when  T=t,  Effectively,  the  time  variable  t is  lost  and  is  not  directly 
recoverable  when  only  the  transient  quadratic  resistance  component  is 
known.  In  none  of  the  correlation  approaches  attempted  was  it  possible 
to  achieve  the  equivalent  of  the  situation  in  Eq.  (47)  where  the  argu- 
ments in  the  two  input  transients  do  not  have  a common  term. 


49 


R- 1928 


What  was  sought  in  attempting  correlation  methods  was  some  means 
of  developing  from  the  observables  a function  of  time  or  time  lags  which 
is  a double  Fourier  transform  of  the  quadratic  frequency  response  func- 
tion. Failure  to  find  such  a function  in  the  present  work  does  not  prove 
that  it  is  impossible  to  do  so.  However  the  nature  of  the  quadratic 
model  is  such  that  it  seems  possible  that  insufficient  information  is 
contained  in  a single  resistance  transient  to  enable  a more  or  less  direct 
identification  of  the  complete  quadratic  frequency  response  function. 
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A FILTERING  APPROACH  TO  SEPARATION 
OF  LINEAR  AND  QUADRATIC  COMPONENTS 

• 

The  unsolved  problems  of  the  last  section  may  be  considered  in 
two  parts:  I)  The  quadratic  component  of  the  resistance  was  not  isolated 

and,  2)  the  quadratic  frequency  response  function  appears  in  an  integral. 

Of  these  two  problems,  the  first  is  at  least  as  critical  as  the  second 
since  the  linear  component  of  the  transient  would  be  expected  to  be  domi- 
nant in  most  experimental  cases.  Even  if  the  integral  can  be  inverted  in 
some  way,  the  linear  component  amounts  to  a "noise"  which  is  large  rela- 
tive to  the  quadratic  signal. 

The  prospects  for  something  useful  resulting  from  a filtering 
operation  on  the  resistance  transient  are  best  discussed  in  conjunction 
with  the  Fourier  Transform,  Eq.  (46).  Figure  14  shows  the  computed  moduli 
of  the  Fourier  Transforms  of  the  various  components  of  the  simulated  results 
for  wave  pulse  No.  1 (Figure  10).  At  the  top  is  the  modulus  of  the  wave 
pulse  transform,  ] H(id)  | . In  the  middle  the  modulus  of  the  linear  compo- 
nent, | G1(u))H(aj)  | , is  shown  as  a solid  line.  Superposed  in  dashes  is 
the  modulus  of  the  transform  of  the  quadratic  component  (the  integral  in 
Eq.46).  Finally  the  modulus  of  the  transform  of  the  total  resistance, 

| R(ui)  | , is  shown  at  the  bottom. 

It  is  clear  from  the  example  in  Figure  14,  as  well  as  the  form  of 
the  integral  in  Eq.  (46),  that  the  quadratic  component  of  resistance  has 
some  frequency  components  the  same  as  those  produced  by  the  linear  part 
of  the  system.  If  all  the  linear  components  are  eliminated  by  some  sort 
of  filtering  a great  deal  of  whatever  information  is  contained  in  the 
quadratic  component  will  also  be  eliminated. 

To  be  more  specific,  suppose  the  resistance  transient  corresponding 
to  Figure  14  is  passed  through  a low-pass  filter  with  cutoff  frequency 
adjusted  so  that  all  frequency  components  above  a non-dimensional  frequency 
of  0.5  are  attenuated  to  much  less  than  1%  of  their  original  magnitude. 

Under  these  circumstances  practically  all  of  the  linear  contribution  to 
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FIGURE  14  MODULI  I OF  COMPUTED  FOURIER  TRANSFORMS 
OF  SIMULATED  RESISTANCE  COMPONENTS, 
WAVE  PULSE  NO.  I 
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the  transform.  Figure  14,  disappears,  as  well  as  the  quadratic  contribu- 
tion to  the  transform  above  a frequency  of  0.5.  What  is  left  of  the 
resistance  transform  will  be  at  frequencies  below  the  lowest  wave  pulse 
frequency,  and  what  will  be  left  in  the  time  domain  will  be  only  a rela- 
tively long  period  excursion  of  resistance.  In  terms  of  the  quadratic 
frequency  response  function  what  is  thrown  away  by  this  operation  is 
defined  by  the  filtering  equation,  Eq.  (28).  Referring  to  Figures  8 and  9 
the  low-pass  filtering  operation  just  described  produces  an  apparent  quad- 
ratic frequency  response  function  which  is  zero  or  negligible  for  | Q2  | 
greater  than  0.5,  and  which  is  appreciable  only  near  the  position  of  the 
mean  added  resistance  operator  (the  axis,  Figs. 8, 9).  So  long  as  the 
low-pass  filter  has  unity  DC  gain,  the  apparent  quadratic  frequency 
response  function  after  filtering  contains  the  mean  added  resistance  opera- 
tor undistorted. 

Thus  a straightforward  low-pass  filtering  operation  appeared  to 
have  some  promise  if  estimates  of  the  quadratic  frequency  response  func- 
tion in  the  neighborhood  of  the  mean  added  resistance  operator  are  all 
that  is  desired.  Because  the  mean  added  resistance  operator  is  the  only 
part  of  the  quadratic  response  function  which  has  found  practical  applica- 
tion, the  approach  appeared  worthwhile  pursuing. 

To  formalize  the  low  pass  filtering  approach  somewhat,  it  is  first 
required  that  a low-pass  filter  frequency  response  function,  Le(uj),  be 
specified  so  that 

G1  (uo)  H(u>)  Le(uu) 

is  zero  or  approaches  zero  for  all  frequencies.  The  filter  should  have 
unity  DC  gain.  Because  the  linear  component  of  resistance  can  have 
appreciable  frequency  components  only  at  frequencies  at  which  the  wave 
pulse  has  appreciable  frequency  components,  the  specification  of  the  cut- 
off frequency  of  the  filter  may  be  made  on  the  basis  of  the  frequency  of 
the  lowest  appreciable  frequency  component  of  the  wave  pulse.  Knowing 
this  frequency,  the  filter  cutoff  must  be  adjusted  so  as  to  attenuate 
resistance  components  at  this  frequency  to  some  very  small  fraction  of 
the  original. 
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Supposing  the  filter  to  be  selected  in  this  way,  an  expression 
for  the  filtered  resistance,  r (t) , may  be  written  with  the  aid  of 
Eqs.  (28),  (45)  and  (46): 

r.(t>  “9oLe(0) 

+ IK  Le (“>) H (^)  H (^)  Exp[  i U)t]dpd0)  (49) 

where  as  in  previous  equations,  corresponds  to  the  sum  or  output  fre- 
quency and  g to  the  difference  frequency. 

The  elimination  of  the  calm  water  resistance,  gQ,  from  Eq.  (49) 
may  be  made  as  has  been  previously  suggested.  This  is  to  fit  a straight 
line  through  beginning  and  end  of  the  filtered  transient  and  subtract  the 
result  from  the  transient  itself.  Then  after  filtering  and  correction 
for  the  calm  water  resistance  what  is  left  may  be  called  r2  (t)  and  is 
written: 

( ^ ) “ JJG2  (^>i42^)Le(a))H(-ii^)H(^|^)  Exp[  iu>t]dpda>  (50) 

In  Eq.  (50)  l-e(u))  is  specified  by  the  filter,  H(uu)  is  the  Fourier  Trans- 
form of  the  encountered  wave  pulse,  and  rg  (t)  is  the  result  of  filtering 
the  observed  resistance  transient  and  correcting  this  result  for  calm 
water  resistance.  Only  the  quadratic  frequency  response  function, 

G2(uUj,<ds)  is  considered  unknown. 

Temporarily  leaving  aside  the  question  of  whether  anything  can  be 
done  with  Eq.  (50),  some  practical  filtering  considerations  may  be  dealt 
with.  The  foremost  question  is  whether  the  specified  filter  can  be 
realized.  In  the  hypothetical  experiment  the  wave  pulse  is  supposed  to 
be  under  sufficient  control  so  that  the  influence  of  frequency  components 
below  some  lowest  encounter  frequency  is  nil.  The  lowest  encounter  fre- 
quency is  determined  by  the  longest  wave  length  assumed  to  be  important 
and  by  model  speed.  In  the  present  example  if  the  longest  wave  length  of 
importance  is  fixed  as  shown  in  Figure  2 and  the  Froude  number  is  varied, 
the  lowest  wave  encounter  frequency  of  importance  varies  only  by  ± 10% 
from  that  shown  in  Figure  14  for  model  speeds  ranging  from  zero  to  twice 
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that  of  the  example.  Accordingly  the  filtering  situation  as  defined  by 
Figure  14  may  be  taken  as  typical.  In  this  case  in  order  to  substantially 
eliminate  frequency  components  above  a non-dimensional  frequency  of  0.5, 
something  like  a six  pole  filter  with  nominal  cutoff  frequency  of  0.25  is 
requ i red . 


If  the  basic  data  is  digitized  before  filtering,  two  methods  are 
available  with  which  filters  of  the  above  specification  may  be  realized. 

The  easiest  is  perhaps  the  FFT  based  fast  convolution.  In  the  present 
example  a recursive  digital  filter  is  difficult  because  of  the  very  low 
cutoff  in  relation  to  the  folding  frequency  range  (in  the  example  the  At 
chosen  results  in  a folding  frequency  30  times  the  required  filter  cutoff). 
However  it  is  feasible  to  implement  a two  stage  digital  filtering  procedure 
where  the  first  stage  removes  enough  high  frequency  response  so  that  the 
data  may  be  decimated,  and  the  second  stage,  operating  on  the  decimated 
data,  achieves  the  desired  final  filter  cutoff. 


Alternately,  there  is  the  possibility  of  incorporating  a real 
analog  filter  into  the  experimental  instrumentation.  The  non-dimensional 
encounter  frequencies  noted  in  Figure  14  are  almost  numerically  equal  to 
real  frequencies  (in  Herz)  for  the  5 foot  model  size  described.  Six-pole 
low-pass  filters  with  fixed  cutoff  frequency  of  0.25  Hz  are  not  now  off- 
the-shelf  items,  but  solid-state  electronic  modules  with  which  such  a 
cutoff  may  be  real ized  are  readily  available.  In  situations  where  the 
model  is  "large",  say  20  feet,  the  filter  cutoff  would  have  to  be  in  the 
neighborhood  of  O.i  Hz.  These  cutoff  frequencies  are  also  considered 
within  reason  for  modern  equipment  and  the  alternative  of  i ncorporat i ng 
real  filters  into  the  experiment  is  also  assumed  to  be  feasible.  In 
such  a case  the  data  digitized  for  analysis  would  be  the  filtered  output. 


Ideally  there  would  be  no  preference  between  digital  and  analog 
filtering  — only  a matter  of  whether  high  frequency  content  would  be 
thrown  away  at  the  outset  or  after  the  experiment  is  over.  In  practice, 
preference  might  be  given  to  analog  filtering  during  the  experiment. 
Judging  by  the  examples  previously  shown,  the  filtered  resistance  excur- 
sion is  apt  to  be  small  relative  to  the  maximum  peak  to  peak  range  of 
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the  total  resistance  transient.  Real  time  analog  filtering  would  make 
possible  the  adjustment  of  gains  so  that  the  best  use  is  made  of  the 
resolution  available  on  the  recording  medium  — a particularly  important 
consideration  if  the  recording  medium  is  analog  magnetic  tape. 


f 

AN  APPROXIMATE  APPROACH  TO  THE  IDENTIFICATION  OF 
THE  MEAN  ADDED  RESISTANCE  OPERATOR 

Presuming  that  the  linear  component  of  the  resistance  transient 
may  be  removed  one  way  or  another  by  filtering  there  is  left  Eq.  (50), 
which  has  the  form  of  a pure  quadratic  system.  The  function  to  be 
identified  is  (^,  - -|) . Continuing  the  approach  in  the  last  section, 
if  the  pass  band  of  the  filter  Lg(uo)  is  reduced  even  further  than  is 
required  to  eliminate  the  linear  component,  the  effective  quadratic  fre- 
quency response  function  is  concentrated  along  the  line  u^=-iu3  (n?=0, 

F i gures  8,9) . 

Considering  the  variation  of  the  actual  quadratic  frequency  response 
function  near  the  line  <^=-0)  , according  to  the  general  symmetry  proper- 
ties of  the  function  it  is  expected  that 


a_ 


u>=0 


[0  + i F ( (p) ] 


in  which  F ( (p)  is  an  even  function  in  p.  Accordingly,  for  values  of  uj 
not  far  from  zero  it  would  be  expected  that  the  quadratic  frequency 
response  function  would  have  the  form: 


Fr(P)  + i U)F ( (p) 

in  which  FR(p)  is  the  mean  added  resistance  operator  [Ga(-|  , - ^)],  and 
is  even  in  p. 

Now  if  it  is  assumed  that  the  low  pass  filter  can  be  made  suf- 
ficiently narrow,  the  effective  quadratic  frequency  response  function 
may  be  approximated  by: 
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K2(^,^)  = La(uj)[FR(P)  + i u>F  j (P)  3 (51) 

In  the  approximation  it  is  assumed  that  L (0)=1  so  that  F_({3),  the  func- 
tion  of  interest,  is  not  distorted.  It  is  further  assumed  that  L (u>)  is 

3 

nearly  zero  for  values  of  uj  for  which  [F  (p)  + iu)F  ((3)]  is  not  a good 

K I 

representation  of  the  actual  frequency  response  function. 

The  filter  in  the  above  approximation  is  assumed  to  be  an  ana  lysis 
filter  which  has  a lower  cutoff  frequency  than  any  real  filter  which  might 
have  been  applied  before  digitizing  the  data.  There  is  no  reason  why  this 
filter  has  to  be  realizable  in  the  real  time  sense,  and  some  convenience 
in  taking  it  as  non-real izable.  Thus  L (m)  will  be  taken  as  real  and 

3 

even  in  uj. 

In  general,  the  approximation,  Eq.  (51)  must  correspond  to  an  effec- 
tive quadratic  impulse  response,  k (tj,^).  Applying  Eq.  (20): 

k8(ti»ta)  = JjExpnVi  + iVa^K,^ 


1 1 

2 (2*)* 


JjExp[i(U  + 


t. -t  UJ+P 

ip  >—' >dujdf3 


(after  making  the  usual  transformation;  0)=^+^,  and  (3=0^-^) 


(52) 


Lett i ng : 


(53) 


and  substititing  Eq.  (51)  into  Eq.  (52),  and  noting  that  F.  (p)  and  F.  (p) 

K I 

are  even  in  P,  and  that  la(uu)  is  even  in  uj,  there  is  finally  obtained: 


mt.> 

+ { vt*>-  mt»> 


(54) 
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where: 


w 

1 

= 2n 

JLg  (u>)  cosT2uxiuj 

(55) 

w 

-1 

2rr 

j ooL a ( oj)  s inTsuxluJ 

(56) 

VV 

1 

= 2n 

JfrO)  cost,  pdf} 

(57) 

MV 

1 

= 2n 

JF,  (p)  cosTj  pdp 

(58) 

The  apparent 

quadi 

-atic  impulse  response  Eq. 

(5*+)  has  the  correct 

symmetry.  Interchange  of  ta  and  t leaves  To  unchanged,  and  Tx  negative, 
but  the  function  is  the  same  because  both  Eq.  (57)  and  (58)  are  even  in  Tj  . 
Equations  (55)  and  (56)  are  just  Fourier  sine  and  cosine  transforms  of  a 
filter  frequency  response  which  is  specified.  Given  what  is  known  about 
the  quadratic  frequency  response  function  the  left  hand  sides  of  Eq.  (57) 
and  (58)  are  expected  to  be  absolutely  integrable  so  that  the  inverses 
may  be  wr i tten: 

Fr(P)  = Jfj  (T^ExpT-ipTjdTj 

= Jf.dJcosT^d^  (59) 

(since  fx (T  ) is  even  in  T ) 
s imi larly: 

F,(p)  = Jf8 (T x ) cosTj pdTt  (60) 

Now  assuming  that  the  analysis  filter  has  been  applied  to  the 
resistance  transient  and  the  calm  water  resistance  has  been  corrected  for, 
the  result  (denoted  r ( t ) ) may  be  written: 

r3  (t)  = JJka(t1,ts)h(t-t1)h(t-ts)dt1dts  (61) 

Equation  (6 1 ) is  exactly  analogous  to  the  third  term  of  Eq.  (18).  Appli- 
cation of  the  Parseval  Theorem  would  result  in  a form  exactly  the  same 
as  Eq.  (50)  with  the  substitution  of  Eq.  (51)  for  the  effective  quadratic 
frequency  response  function. 
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Substituting  Eq.  (5*0  in  Eq . (6l)  and  making  the  variable  trans- 
formation defined  by  Eq.  (53)  there  results: 

<3  (t)  - Jf4  (T1)j£a(T2)h(t-Ts-T1)h(t-Ts+Ti)dT2dT1 

+ Jf8(Ti)K(T8)h(t"T8“Ti)h(t-|?+TiJdTadTi  (62) 

Alternately,  by  replacing  the  quadratic  frequency  response  in  the  third 
term  of  Eq.  (46)  by  the  approximation,  Eq . (51)  there  results  the  Fourier 
transform  of  r3  (t) : 

Rn  (u))  = Jl-a(uO  [FR(P)+ iu)F(  (6)  ]H(^|^)H(ji^)d3  (63) 

In  either  form  the  result  of  the  approximation  to  the  apparent 
quadratic  frequency  response  is  to  separate  an  unknown  function  of  two 
variables  into  the  product  of  a known  function  and  an  unknown  function  of 
one  variable.  Equations  (62)  and  (63)  are  equivalent  in  the  sense  that 
application  of  the  Parseval  Theorem,  Eq.  (35),  to  Eq.  (62)  yields  Eq.  (63). 

A start  at  a practical  data  reduction  method  based  upon  Eq.  (62) 
or  (63)  was  first  made  by  considering  Eq.  (62),  the  time  domain  version. 
Because  the  simulated  transients  correspond  to  response  to  a wave  eleva- 
tion forward  of  the  model  it  is  known  from  the  work  of  Refs.  1 and  6 
that  the  imaginary  part  of  the  function  near  the  mean  added  resistance 
operator  will  vary  strongly,  that  is,  F ( (6)  will  be  relatively  large. 
However,  it  was  shown  in  Refs.  1 and  6 that  a simple  time  shift  of  input 
relative  to  output  tends  to  minimize  the  magnitude  of  the  imaginary  part 
of  the  function  near  the  line  0^=-^.  In  the  case  of  the  present  simula- 
tion the  actual  best  time  shift  is  the  same  as  that  in  Ref.  6,  or 
18  points.  In  order  to  simplify  a first  trial  it  seemed  reasonable  to 
assume  that  the  shifting  operation  would  make  F ( (p)  very  small  (that  is, 
the  imaginary  part  of  the  effective  quadratic  frequency  response  function 
would  be  nil).  Under  this  further  approximation  fa  (Ta ) , (Eq.58)  is  also 
negligible  and  the  filtered,  time  shifted  resistance  response  was  approx- 
imated as: 

r4(t)  * Jfi(Ti)Iia<T2)h(t-Ta-T1)h(t-T8+T1)dTadT1  (64) 
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For  purposes  of  an  actual  digital  computation  it  was  assumed  that 
the  same  type  of  representation  which  had  produced  the  simulated  trans- 
ients (Eq.29)  would  be  adequate.  Thus  to  transpose  to  the  digital  form: 

t = pAt 
Tj  = jAt 
T = kAt 

ond  Eq.  (64)  becomes: 


r (p)  = E S f , . h(p-k-j)h(p-k+j) 


j k 


lj  ak 


(65) 


where  h(p-k-j)  - h(p-k+j)  represent  points  on  the  time  history  at  time 
steps  (p-k-j)  and  (-— k+j) . Similarly,  r4(p)  is  the  resistance  time  history 
at  time  step  p„  In  Eq.  (65)  discrete  approximations  to  f (T, ) and  i (T  ) 
were  made  as  follows: 

MTi)  = f ]j  6(Ti-jAt)  (66) 

£a(Ts)  = *aj  6(Ts-kAt)  (67) 

Now  taking  the  Fourier  transforms  of  Eq.  (66)  and  (67);  there  results: 

FrO)  = f ,0  + 2 I f j j cos  (jfJAt)  (68) 

j 

La((«)  = iao  + 2 Z ia|<cos(kuAt)  (69) 

k k-1,2... 

Manipulating  Eq.  (65)  slightly, 

r4(P)  = E f,  . C„.  (70) 


: >J  PJ 


j =0, 1 , 2. . . 


where 


C„i  = 4:  1 ^akh(P“k-j)h(p-k+j) 


PJ 


■*  k=-n 


(71) 


with  q = 1 


q.  = 2 for  j ^ 0 
J 
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Once  the  filter  impulse  response  {I  ,)  is  specified,  Eq.  (71)  is  calcul- 
able  as  a function  of  j for  each  response  time  step,  p.  Accordingly, 

Eq.  (70)  represents  a series  of  linear  algebraic  equations  in  the  unknown 
discrete  impulse  response,  fjj'  It  was  hoped  that  this  series  of  equa- 
tions could  be  solved  in  a least  squares  sense  for  the  f^,  and  thus 
ultimately  estimates  could  be  made  of  F^(p)  via  Eq.  (68). 

The  steps  done  in  trying  this  approach  out  were  as  follows: 

a.  Compensate  the  simulated  resistance  transient  for 
calm  water  resistance. 


b.  Assume  a Hanning  Type  Filter  and  compute  1 


ak' 


= — [ 1+cosTTk/m] 


(72) 


where: 


£ [l+cosnk/m] 
k=-m 


and  m is  chosen  to  control  the  cutoff. 

Specifically,  the  corresponding  frequency 

response,  L (uj)  is  unity  for  <u=0  and  roughly 
3 

1/2  for  a*=TT/nv^t. 

c.  Compute  L (oj)  , Eq.  (69). 

3 

d.  Filter  the  resistance  transient  with  L (uj) 

3 

utilizing  the  FFT  fast  convolution  method. 

e.  Shift  the  wave  pulse  time  series  relative  to 
the  resistance  time  series  as  previously  noted. 

f.  Compute  C . for  each  point  in  the  wave  time  series, 
and  accumulate  the  coefficients  in  a least  squares 
fit  of  Eq.  (70) . 

g.  Scale  the  coefficients  in  the  resulting  L.S.  fit 
equations  as  required. 

h.  Solve  for  the  f^  using  a standard  Gauss  elimination 
technique  (Subroutine  GELG  or  DGELG,  IBM  Scientific 
Subroutine  Package) . 
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i.  Estimate  F^(fi)  via  Eq.  (68)  using  the  results  of 
step  h. 

The  above  procedure  was  applied  to  the  simulated  sets  of  transient 

data  with  various  choices  of  m,  that  is,  various  choices  of  low  pass 

filter  cutoff.  The  number  of  estimates  of  f,.  solved  for  was  varied  from 

lj 

10  to  30.  In  no  case  were  the  results  of  the  procedure  meaningful. 

Double  precision  accumulation  of  L.S.  fit  coefficients  and  solution  for 
fjj  was  necessary  in  order  to  get  any  sort  of  solution  for  the  case  that 
the  scaling  (step  g)  involved  only  division  of  the  entire  set  of  equations 
by  the  largest  absolute  value  in  the  coefficient  matrix.  A scaling  pro- 
cedure was  devised  with  which  all  the  coefficients  on  the  left  hand 
coefficient  matrix  were  of  the  order  of  unity  and  by  which  each  row  and 
column  of  the  coefficient  matrix  contain  a value  near  unity.  The  solution 
after  this  scaling  was  no  more  meaningful  than  those  mentioned  previously. 

The  coefficient  matrix  obtained  in  step  f of  the  procedure  is 
certainly  numerically  ill-conditioned.  Those  answers  which  were  returned 
for  the  f^  were  typically  alternately  positive  and  negative  and  of  about 
the  same  magnitude  irrespective  of  the  magnitude  of  j . To  see  what  answer 
should  have  been  obtained,  the  simulated  added  resistance  operator,  Fig- 
ure 8,  was  transformed  numerically  in  accordance  with  Eq.  (57)  and 
evaluated  at  a time  step  corresponding  to  that  used  in  the  analysis. 

These  results  indicated  that  the  absolute  values  of  the  f,.  should  have 

'J 

decreased  as  j increased,  the  magnitude  should  have  been  small  enough  to 
disregard  after  30  time  steps,  and  that  semi-oscillations  plus  to  minus 
should  have  been  about  4 time  steps.  This  evidence  along  with  the  nature 
of  the  fitted  results  suggested  an  over-fit.  Thus  i 1 1 -cond i t i on i ng  appears 
not  to  be  merely  numerical  but  that  there  is  really  not  enough  independent 
information  left  after  the  heavy  filtering  operation  to  enable  a solution 
for  the  30  or  so  values  of  f^.  which  seem  to  be  required  for  a reasonable 
estimate  of  the  mean  added  resistance  operator. 

There  is  of  course  the  possibility  that  the  assumption  that  Fj((3) 
is  negligible  is  at  fault.  However,  the  low  pass  filtering  operation 
plus  the  shifting,  was  certainly  enough  to  make  the  imaginary  part  of 
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the  apparent  quadratic  frequency  response  function  relatively  small, 
if  not  negligible.  Given  that  the  results  obtained  under  the  assumption 
of  negligible  F (p)  were  not  even  in  the  ballpark  it  was  doubted  that  a 
procedure  involving  Eq.  (62)  rather  than  Eq.  (64)  would  be  worth  the 
substantially  increased  effort,  and  accordingly  the  above  approach  was 
abandoned . 

Some  effort  was  made  to  develop  a similar  procedure  with  Eq.  (63), 
the  frequency  domain  version  of  Eq.  (62).  Referring  to  Eq.  (63),  the 
product : 

L>)H(^)H(^) 


is  calculable  from  the  Fourier  transforms  of  the  wave  pulse.  In  addition 
to  a relative  time  shift,  a re-location  of  the  (arbitrary)  position  of 
time  = 0 was  advisable  before  carrying  out  the  Fourier  transform  on  both 
the  wave  and  resistance  transients  so  that  the  arguments  of  each  were 
slowly  varying.  The  net  result  of  the  development  was  an  estimation 
equation  of  the  form: 

RgU)  = JfrO)  Cx  (u>,p) 


+ Jf,(P)  CsK(3) 

+ 1 Jfr(p)  CsKp) 

+ i Jf ( (P)  C4(u>,p)  (73) 


in  which  the  Cn(u>,P)  are  calculable  functions  of  m and  p.  This  approach 
was  carried  through  in  a similar  manner  to  that  indicated  in  Eq . (65) 
through  (71).  The  results  obtained  after  applying  the  approach  were  just 
as  bad  as  before  --  regardless  of  whether  F j (p)  was  assumed  to  be  neglig- 
ible or  not.  The  basic  problem  with  this  approach  is  that  after  applying 
the  Fast  Fourier  Transform  to  the  resistance  transient  and  carrying  out 
the  low-pass  filtering  there  are  a limited  number  of  values  of  R^  (uo) 
available.  In  the  present  case  with  IK  arrays  defining  the  transient, 
and  a filter  narrow  enough  to  eliminate  the  linear  component  of  resist- 
ance, there  were  only  a dozen  or  so  estimates  of  R3  (uj)  which  were  not 
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negligible  due  to  the  effect  of  the  low-pass  filter.  Thus  FD(ji)  and  F (3) 
needed  to  be  approximated  by  no  more  than  a dozen  discrete  values  along 
the  £5  axis  from  zero  the  the  folding  frequency.  This  frequency  resolu- 
tion is  much  too  coarse  for  a practical  data  reduction  procedure.  On 
paper,  the  solution  would  be  to  add  zeros  to  the  transient  arrays, 
doubling  their  length  so  that  more  estimates  could  be  obtained  in  the 
pass-band  of  the  filter.  This  approach  is  thought  illusory  and  was  not 
attempted  because  all  that  is  achieved  is  an  interpolation  — no  new  and 
independent  information  is  achieved.  Just  as  in  the  time  domain  approach, 
the  few  results  achieved  with  the  frequency  domain  approach  indicated  that 
there  is  really  not  enough  independent  information  left  after  filtering  to 
allow  reasonable  answers. 
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CONCLUDING  REMARKS 

Because  the  experimental  determination  of  the  characteristics  of 
mean  ship  model  resistance  added  by  waves  is  one  of  the  more  difficult 
and  time  consuming  problems  in  seakeeping  towing  tank  practice,  it  was  of 
practical  interest  to  see  if  a technique  involving  wave  pulses  could  be 
developed,  and  if  so,  if  it  promised  better  efficiency. 

Toward  this  end  a fairly  realistic  digital  simulation  was  made  of 
the  head  wave  pulse  experiment  and  of  the  resulting  resistance  transient. 
There  appears  to  be  little  doubt  that  such  an  experiment  is  feasible,  and 
that  the  influence  of  calm  water  resistance  can  be  eliminated  in  a straight- 
forward way  because  the  "memory"  of  the  nonlinear  part  of  the  added 
resistance  is  not  exorbitantly  long. 

The  development  of  a data  reduction  procedure  was  approached  from 
the  point  of  view  that  any  practically  attractive  data  reduction  procedure 
must  involve  a very  few  wave  pulse  runs,  preferably  only  one.  This  point 
of  view  was  re-inforced  by  the  wave  pulse  simulation  because  it  appeared 
that  the  detailed  shape  of  the  wave  pulse  would  be  seriously  affected  by 
small  changes  in  wavemaker  and  model  control  parameters.  For  example, 
wave  pulse  time  histories  of  exactly  the  same  shape  but  different  ampli- 
tudes appear  to  require  more  precise  control  of  experiment  than  is 
ordinarily  possible.  Accordingly  it  appeared  that  if  reasonable  estimates 
could  not  be  obtained  from  one  experiment,  there  would  be  little  point  in 
attempting  methods  which  required  multiple  experiments  for  any  purpose 
other  than  conf i rmat ion, or  for  improving  quality  of  estimates  through 
simple  averaging. 

Several  correlation  and  Fourier  analysis  approaches  to  the  ident- 
ification from  transient  data  of  the  quadratic  frequency  response  function 
for  added  resistance  were  investigated.  All  the  approaches  were  analogous 
to  one  or  another  of  the  known  approaches  for  outputs  of  linear  and  quad- 
ratic systems  in  response  to  random  excitation.  None  of  these  investigations 
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resulted  in  a clean  approach  to  the  problem.  Evidence  was  developed  that 
a conceptually  clean  approach  may  be  impossible.  With  respect  to  a 
recovery  of  the  ent'i  re  quadratic  frequency  response  function  for  added 
resistance,  it  appears  that  insufficient  information  is  contained  in  a 
single  transient.  It  should  be  remarked  that  this  is  in  contrast  to  the 
case  for  random  Gaussian  excitation.  In  this  latter  case  a clean  approach 
arises  and  is  traceable  to  statistical  properties  of  the  expected  values 
of  products  of  input. 

By  taking  advantage  of  the  general  properties  of  the  quadratic  model 
which  appears  to  be  valid  for  added  resistance,  it  was  possible  to  develop 
a "dirty"  approach  for  the  identification  of  the  mean  added  resistance 
operator  from  a single  transient  experiment.  However  no  meaningful  results 
could  be  obtained  with  this  latter  approach.  Those  results  which  were 
obtained  imply  that  there  is  not  enough  information  in  a single  transient 
experiment  to  enable  the  identification  of  even  a special  portion  of  the 
quadratic  frequency  response  function  (the  mean  added  resistance  operator). 

It  appears  that  if  there  is  a practical  wave  pulse  technique  for 
f added  resistance  experiments,  it  is  not  yet  in  hand. 
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PRINCIPAL  NOTATION 
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L1(a>),  L2(uj) 
La(oj),  Le(ou) 

L, 
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W<  VT.> 


approximating  functions 

impulse  responses 

weighting  coefficients 

linear  frequency  response  function 

quadratic  frequency  response  function 

gravitational  constant 

zeroth  degree  kernel  (calm  water  resistance) 

linear  kernel  (surge  force  impulsive  response) 

quadratic  kernel  (added  resistance  impulse  response) 

Fourier  transform  of  wave  pulse 

encountered  wave  pulse 

wave  pulse  at  tank  position  x 

i nd ices 

linear  frequency  responses  (filters) 

linear  weighting  coefficient 
impulse  responses 
quadratic  weighting  coefficients 
resistance  transient 
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t ime  var iables 
t ime 

model  velocity 
tank  pos  i t ion 
difference  frequency 
time  interval 
dummy  time  variables 
difference  frequency 

sum  frequency  (equals  encounter  frequency) 
encounter  frequency 
wave  frequency 

frequency  of  wave  of  ship  length 
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